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Q ■ Abstract 

o : 

(N ; 

^ ! Finding the eigenstates of the total Hamiltonian H or its diagonalization is the impor- 

^ ■ tant problem of quantum physics. However, in relativistic quantum field theory (RQFT) 

' its complete and exact solution is possible for a few simple models only. Unitary trans- 
formations (UT's) considered in this survey do not diagonalize if, but convert H into a 

^ ' form which enables us to find approximately some H eigenstates. During the last years 

^ I there have appeared many papers devoted to physical applications of such UT's. Our 

^ . aim is to present a systematic and self-sufficient exposition of the UT method. The two 

5 ■ general kinds of UT's are pointed out, distinct variations of each kind being possible. We 

O I consider in detail the problem of finding the simplest H eigenstates for interacting mesons 

^— > . and nucleons using the so-called "clothing" UT and Okubo's UT. These UT's allow us to 

Zj ■ suggest definite approaches to the problem of two-particle (deuteron-like) bound states 

pH I in RQFT. The approaches are shown to yield the same two-nucleon quasipotentials in 

>• . the first nonvanishing approximation. We demonstrate how the particle mass renormal- 

^ ' ization can be fulfilled in the framework of the "clothing" procedure. Besides the UT of 

H : 

, the Hamiltonian we discuss the accompanying UT of the Lorentz boost generators. 

1. Introduction 

The so-called unitary transformation (UT) method has the same age as the quantum 
theory itself. Its first applications for constructing Hermitian effective interactions (HEI) 
can be found in and 0] (in this connection, see a review article 0] where the differ- 
ent perturbation expansions of HEI were discussed). In quantum field theory the first 
considerations using the UT method were given by WentzelQ and Heitler|^. 

A number of the schemes for reduction of the exact eigenvalue problem to the model- 
space problem via the various UT's were put forward in 50's (0,[0,0) during the exten- 
sive development of the meson theory of nuclear forces (see,e.g.,0). Owing to the work 
||T0| this approach has proved to be very useful in studies of electromagnetic (e.m.) inter- 



actions with nuclei for nonmesonic channels (in particular, when constructing the effective 



1 



operators of meson exchange currents /see, e.g., |TT]], surveys [O, [O and refs. therein/). 

Along with this guidehne the UT method was used for the formulations of RQFT in 
terms of physical or "clothed" particles [|l^, [|16], Jl^, ||18|, see also our talks at the 

recent conferences |T9| and pO|). 

By using the Okubo UT the authors of |2T] constructed effective generators for the 
Poincare algebra, acting on nucleonic degrees of freedom only. It was made in the frame- 
work of perturbation theory for a simple model of "spinless" nucleons exchanging scalar 
mesons. Therefore, we have an instructive example of how some noncommuting Hermitian 
operators can be reduced by one and the same UT to Okubo's block form. 

Late 80's and 90's have brought a renewed interest to this area. First of all, we mean 
applications ([^,p3|) of the UT method to a covariant treatment of the two-body bound- 
state problem (cf. p^). The corresponding transformed Hamiltonian and boost operators 
do not couple (in the second order in meson-nucleon coupling constants) the nucleon 
(no-meson) subspace with its complement in the full Fock space of hadron states. In the 



work p5[ the same method has been employed to derive the effective nucleon-nucleon and 
nucleon-antinucleon interactions starting from a field Hamiltonian with the exchange of 
TT, p, uo and a mesons. Then, within the Hartree approximation, these effective interactions 
have been introduced for describing the saturation properties of nuclear matter. Note also 
recent explorations ^ of the hard-core problem in the theory of nuclear forces. 

29| ) some extensions of the UT method have been suggested for con- 
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structing effective current operators in the theory of photomeson processes on nuclei (see 
also ref. where one can find the calculations of the four structure functions for pion 
electroproduction on the deuteron near threshold). Recently,the UT method has been 
used [0 for deriving effective two-particle one-meson exchange potentials in the instant 
and front forms of relativistic quantum mechanics. At last, with the aid of a modification 
of the method the authors of have proposed a meson-exchange model for vrN scat- 
tering and 7N-7rN reaction. Certainly, one may say that nowadays the UT method has 
survived its second birth. 

This review is focused upon an application of UT's in RQFT and aimed at an approx- 
imate treatment of the physical vacuum, the observable one-particle and two-particle 
bound and scattering states. These states have a common feature, viz., they do not 
change in time. The simplest example is the state without observable particles, i.e., 
the physical vacuum. Other examples give free particles with a definite four-momentum 
(elementary particles, atoms and nuclei in their ground states). 

Our point of departure is that such states should be eigenvectors of a Hamiltonian if, 
which are stationary. In the context, the "bare" states, i.e., eigenvectors of a free part Hq 
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of the total Hamiltonian H, are inappropriate for theoretical description of the physical 
objects. Firstly, they change in time although the possible transitions to other states are 
virtual and can be considered undetectable because of their small probabilities. Secondly, 
Hq has no eigenvectors which would correspond to bound states (e.g., the hydrogen atom 
or the deuteron). 

From this postulate it follows such a definition for a bound state of the deuteron- 
type in RQFT : it should be described by a proper eigenvector of the total Hamiltonian 
of the theory. In other words, the deuteron problem should be reduced to an exact or 
approximate solution of the H eigenvalue problem. There are various approaches to this 
problem, viz., within the Bethe-Salpeter formalism, via its three-dimensional versions. 



etc. (see, e.g., |^3|, [Q, and refs. therein). However, the links between the respective 
" wavefunctions" (for instance, the Bethe-Salpeter amplitudes for bound states) and the 
exact or approximate H eigenvectors are very sophisticated. According to our postulate 
similar links in the framework of the approach developed below are much more direct and 
transparent (see Sect. 4). In addition, the description of particle scattering can be reduced 
to calculation of the relevant scattering wavefunctions. 

The reasons for employing just RQFT when describing bound and scattering states 
are well known. Being not satisfied with the multitude of disconnected phenomenological 
explanations we strive for a unified description of Nature. RQFT's are the best known 
candidates for unified theories. Firstly, they give a qualitative and natural consideration of 
particle creation and destruction. Secondly, local RQFT's ensure in a sense the relativistic 
causality unlike phenomenological approaches. 

General idea of constructing the UT's in question can be formulated in the following 
way. 

Any UT of H can be considered as a transformation expressing the iY-matrix deter- 
mined with respect to a new basis in the Hilbert (Fock) space of a given physical system 
through the if-matrix with respect to the old basis in the same space (details see in 
Sect. 6). Normally, the latter is composed of eigenvectors (the "bare" states) of the 
free part Hq in the partition H = Hq + V , where V represents the interaction between 
the fields involved. In general, the new if-matrix turns out to be more complicated. 
For example, in the case with Yukawa coupling V ^ J il!{x)'y5tp{x)(f){x)dx, that con- 
serves the baryon number and may change the meson number merely by unit, nonzero 
elements of the initial iZ-matrix are either diagonal, or near diagonal (e.g., they can be 
of the (A^' I H I 7rA^)-kind). In a new representation for this matrix all its elements 
can be nonzero (e.g., along with the aforementioned elements one can meet the elements 
(A^'A^' I H I A^A^)), (vr'A^' | H \ nN), (N' \ H \ ttttN), and the others in which the meson 
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number can be arbitrarily altered^). The elements {N'N' \ H \ NN) contribute to the 
A^A^ scattering even in the first order of perturbation theory. Of course, there are other 
elements of the transformed if-matrix, which contribute to the A^A^ scattering as well, 
but in higher orders. 

Let us impose the following constraint upon the transformation associated with the 
basis change: in the new if-matrix the A^A^ interaction should be described only through 
the elements {N'N' \ H \ NN), i.e., all the other elements which could contribute to the 
A^A^ scattering should be zero. This means that the if-matrix must be reduced to a 
block diagonal form where specific off-diagonal blocks consist of zero elements (other 
clarifications can be found in Sect. 6). This requirement (the Okubo condition0) can be 
replaced by the other constraints, more flexible and easier realizable. 

Our aim is to present a systematic exposition of the UT method. A main attention 
within the method is paid to the two approaches, viz., the clothing procedure and the 
blockdiagonalization after Okubo. We strive for a self-sufficient presentation which may 
be understood without referring to original papers. As a rule, we avoid to point out 
mistakes and obscurities in the latter (but sometimes allow ourselves to note their tacit 
assumptions) . 

This review is organized as follows. 

In Sect. 2 we consider the problem of obtaining the simplest H eigenstates. The lowest 
H eigenstate fl can be juxtaposed to the state without observable particles (the physical 
vacuum). Further, we seek one-particle-like H eigenstates which have "bare" partners 
(e.g., in the case of the interacting pion and nucleon fields the "bare" one-meson a^{k)Qo 
and one-nucleon b'^{p)Qo states, Qq being the "bare" vacuum). These H eigenstates may 
be called "clothed" because meson-nucleon interaction is taken into account when 
constructing such states. At this point, instead of the usual "bare" creation-destruction 
operators a^{k), a{k), ... we introduce the "clothed" operators aj{k), adk), so that 
the physical one-meson state is described by the vector ac^{k)VL. 

While finding the above H eigenstates one has to determine a "clothing" UT such that 
the transformed Hamiltonian does not contain the interaction terms which correspond to 
some virtual energy-nonconserving processes (e.g., A^ — > vrA^, tt NN ). These terms 
are called "bad". 

As one might expect the mass of the clothed particle turns out to be inequal to the 
respective bare mass. The related problem of mass renormalization is explored in Subsect. 
2.5. 

The clothing transformation of the Hamiltonian results in its representation through 
^In order to distinguish the matrix elements with respect to the new basis we employ round brackets. 
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the clothed creation-destruction operators. A basic task of Sect. 3 is to investigate how 
transformed generators of the Lorentz boosts depend on the clothed operators. We show 
that the bad terms can be removed via the same clothing transformation simultaneously 
from the total Hamiltonian and these generators. A consequence of such reduction is that 
the "clothed" vacuum fl remains invariant under Lorentz transformations (unlike the 
bare vacuum Qq) while the clothed one-particle states have the proper transformation 
properties (in particular, their momenta being suitably changed). 

Along with the clothed one-particle states the total Hamiltonian can have one- 
particle-like states (e.g., bound states of the deuteron-type) without any bare partners. In 
Sect. 4 within our clothing procedure we suggest an approximate way of finding such states. 
The resulting bound-state equations resemble the Schroedinger equation for stationary 
states in nonrelativistic quantum mechanics, where the usual potentials are replaced by 
the interactions between clothed particles (the quasipotentials) . Explicit analytical ex- 
pressions for the model nucleon-nucleon and meson-nucleon quasipotentials are given in 
Subsect.4.3. 

In Sect. 5 we discuss some modifications and extensions of the clothing approach. 

Sect. 6 starts with a definition of the similarity transformation H — > U~^HU {U'^ — 
W) as a transformation of the matrix {n' \ H \ n) (determined with respect to a set of 
orthonormal vectors | n)) into the matrix {u' \ H \ u) with respect to another basis | v). 

The Okubo condition on U and the corresponding decoupling equation from [7] are 
presented in Subsect. 6.2. 

Two operator kinds of the Okubo UT are considered in the framework of the Okubo 
approach using the new particle creation-destruction operators which correspond to the 
clothed operators determined in Sect. 2. This allows us to establish the relation to the 
clothing procedure and its modifications (cf. Sects. 2 and 5). 

Some original results of this work are summarized in Sect. 7. 

More technical details and some auxiliary calculations are deferred to Appendix A. 
Appendix B is devoted to some mathematical aspects related to the UT method. We 
suggest an algebraic approach which enables to treat the clothing UT as an element of an 
algebra lacking any operator representation. At last. Appendix C exemplifies an explicit 
solution of Okubo's decoupling equation for a simplified field model. 

2. Clothed Particles in Quantum Field Theory 

The notion of clothed particles will be considered using the following model: a spinor 
(fermion) field -0 interacts with a neutral pseudoscalar meson field (p by means of the 
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Yukawa coupling. The model Hamiltonian is H = Hq + V where 

Ho = J ^Pix)[-i-fV + mo]i:{x)dx + [n^ix) + V(p{xf + fio^cp^^x)] dx (2.1) 

V = ig J tp{x)^5^{x)(j){x)dx. (2.2) 

For simplicity, we do not employ a more refined form of H properly symmetrized in the 
fields involved (see, e.g., JlSl, |3^). This model has much in common with more realistic 
models for the interacting fields (e.g., the nucleon isodoublet (p,n) interacting with the 
meson isotriplet (7r+, 7r°, vr"). 

The Hamiltonian can be expressed in terms of bare destruction (creation) operators 
a{k) (a^(fe)), b{p,r) (6^(p, r)) and d{p,r) {d^{p,r)) of the meson, the fermion and the 
antifermion (see Eqs.(2.8) and (2.16)). Here k and p denote the momenta; r is the spin 
index. An exact definition of the bare operators which we use will be given in Subsect. 
2.2. In what follows, the set of all these operators is denoted by a symbol a, while Op is 
used for one of them. The state without bare particles Qq and the bare one-particle states 
a^{k) Qq, h'^{p^r) VLq and d'^{p,r) VLq are not H eigenvectors. 

2.1 Clothed Particle Operators and States 

Now, we introduce new destruction(creation) operators 

ac{k){ac\k)), hc{p,r){h^\p,r)) and dc{p,r){dj{p,r)) yk,p,r (2.3) 

with the following properties: 

i) The physical vacuum (the H lowest eigenstate) must coincide with a new no-particle 
state fl, i.e., the state that obeys the equations 

ae(fc) I Q) = b,{p, r)\Q) = d,{p, r)\n)=0 Wk, p, r (2.4) 

= 1 

ii) New one-particle states aJ{k)Q etc. are H eigenstates as well 

iii) The spectrum of indices that enumerate the new operators must be the same as 
that for the bare ones (this requirement has been used when writing Eq.(2.3)) 

iv) The new operators satisfy the same commutation rules as do their bare partners. 
For instance, 

K(fc),a,t(fc')] = 5ik-k') , {Kip,r),bJip',r')} = {d,ip,r), dj {p' ,r')} = 5„.'6{k-k') . 

(2.5) 
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Following (Ch.XII) we shall call clothed the new operators and states. Note 

that the name is sometimes used in a sense which differs from that defined by the points 
i) - iv). 

As one can see, the problem of clothing is equivalent to determination of some H 
eigenvectors. In fact, the property iii) means that we do not pretend to find all H eigen- 
states which are one-particle-like. For example, H may have a deuteron-like eigenstate 
with a mass < 2m, where m is the nucleon mass. No bare one-particle state coresponds 
to such a state. Now we intend to find only those one-particle-like eigenstates of H which 
have bare partners. 

One should stress that the clothing problem may turn out to be unsolvable. A solv- 
ability condition will be pointed out in Subsect 2.4. Note also that the properties i) - iv) 
can be supplemented by some physical constraints which will be discussed in Sect. 3, but 
are not needed here. 

Some clothing procedures have been realized within simple field models (see, e.g., 
1^ , |15|, [l^). In the paper, we use a kind of perturbation theory developed in |l^ and 



18| . It can be applied to any field theory to yield an approximate solution of the problem. 

2.2 Bare Particles with Physical Masses 

By definition, the bare one-fermion eigenstate | p,r)o of the operator Hq, being 
simultaneously the eigenstate of total momentum P, belongs to the Hq eigenvalue 
Ep = \^p'^ + mo^. Let us consider an H eigenstate | p,r) for which | p, r)o is a zeroth 
approximation (ZA). Perturbation theory shows that the corresponding H eigenvalue Ep 
differs from Ep. In the relativistic case the function Ep must be of the form -y/p^ + m^ 
where m is the mass of an observed free fermion. We shall call it physical mass. Anal- 
ogously, one can argue appearance of the meson physical mass fi which differs from the 
the trial mass fiQ. 

So, we expect that the physical fermion and meson masses m and fi arise in a natural 
manner when finding H eigenvalues which correspond to the clothed one-particle states. 

Such an introduction of the masses m and fi can be used to divide the total Hamiltonian 
into the new free part Hp and the new interaction Hj. Namely, let us rewrite H = Hq + V 
as if = Hp + Hi where 

Hp = J ij{x)[-ijV + m]^{x)dx + ^ J [n\x) + V^{xf + fi'^(p'^{x)] dx (2.6) 
Hi = V + {mo-m) J i!{x)tlj{x)dx + ^{fiQ^ - /i^) J (j)^{x)dx + M^en , (2.7) 



The decomposition H = Hp + Hi is the well-known trick(see, e.g., [0), but this is 
not necessary for our clothing program: all the following results can be obtained without 
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the introduction of the mass counterterms M^er^- However, it simphfies the program 
reahzation. In other words, the separation H = Hp + Hj may be justified not ab initio 
but post factum. 

The operator Hp can be brought to the "diagonal" formQ 

Hf = J u;}^a^{k)a{k)dk + J EpYP{p. r)b{p, r) + d^{p, r)d{p, r)]dp (2.8) 
by means of the standard expansions 

(f){x) = {2n)-^^^ J {2ujj^)-^^^[a{k) + a^{-k)]exp{ikx)dk (2.9) 

7i{x) = -z(27r)"^/2 J {ujf^/2y^^[a{k) - a\-k)]exp{ikx)dk (2.10) 

ij{x) = (27r)"^/2 f {m/Epf^ I][m(P, r)h{p, r) + v{-p, r)d\-k, r)]ex^{ipx)dp (2.11) 
J J. 

where u{p,r) and v{p,r) are the Dirac spinors, which satisfy the conventional equations 
(^ — ■m)u{p,r) = and (^ + m)v{p,r) = with ^ = Ep'j^ — p^. In these formulae 
Ep = a/p^ + m'^ and uj^ = \Jk^ + /i^. 

2.3 The Unitary Transformation 

The operators (2.3) are the corner-stone of the clothing procedure. Our aim is to find 
clothed operators which should satisfy the requirements i)-iv). Now, the symbol a will 
be used for set (2.3) with being one operator of the set (cf. a and Op). In order to 
implement the properties iii) and iv), we suppose that the clothed operators a are related 
to bare ones a via a unitary transformation 

ap = W^apW, W^W = WW^ = 1 , (2.12) 

where is a function of all the bare operators a. Therefore, Eq. (2.12) represents ap as 
a function(functional) of a. 

Note that W is the same function of either clothed or bare operators (see [l^). Indeed, 
if f{x) is a polynomial or a series of x, the relation /(a) = (a) f {a)W {a) follows from 
Eq.(2.12). Replacing /(a) by W leads to 

W{a) = W^{a)W{a)W{a) = W{a) , (2.13) 

i.e., to the above statement. Hence, the operator ap, when expressed in terms of a, is 
given by 

ap = W{a) ap W\a) . (2.14) 
simple example of the calculation of radiative correction to particle "bare" mass can be found in 
App. C 

■^Nonessential c-number terms are henceforth omitted 
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Unitarity of W is automatically ensured if W is represented as the exponential of an 
antihermitian operator R: W = expR. For a given R, the r.h.s. of Eq.(2.14) can be 
evaluated with the help of 

e^Be-^ = B + [A,B] + ^[A, [A,B]] + ^[A, [A, [A,B]]] + ... (2.15) 

and the commutation rules (2.5). 

In the context, the total Hamiltonian can be written as, H = H{a) = Hp + Hj where 
Hp^a) is determined by Eq.(2.8) and Hi = V{a) + Mren{o,) with Q 

1Q f Tfl 

^(") = (n \V2 / ^P' 1^ T? 77 M/2 ^(P+fe~P0 {^^(pV)75M(pr) fet(pVXpr) + 

(27r)'^/^ J {2uf^EpEp)'^/^ 

+ u{p'r')'-f5v{—pr) h\p'r')d\—pr) + v{—p'r')'-f<^u{pr) d{—p'r')h{pr) + 

+ v{-p'r')-i.,v{-pr) d{-p'r')d\-pr)] [a{k) + o){-k)] . (2.16) 
By using Eq.(2.14), one can replace the bare operators by the clothed ones 

H{a) = H{W{a)aW^{a)) = K{a) . (2.17) 

The operator K{a) represents the same Hamiltonian, but it has another dependence on 
its argument a compared to H{a). K{a) can be found as follows. First, Eq.(2.17) can be 
written as 

K{a) = W{a)H{a)W^{a) . (2.18) 
Second, putting W{a) = expi?(a) and using Eq. (2.15) we have 

H = K{a) = e^ [Hp + Hj] e~^ = 

= Hp{a) + Hj{a) + [R, Hp] + [R, Hj] + ^[i?, [R, Hp]] + ^[i?, [i?. Hi]] + ... (2.19) 

Eq. (2.19) gives a practical recipe for the K{a) calculation: at the beginning one replaces 
a by a in the initial expression H{a) and then calculates W{a)H{a)W\a) using Eqs. 
(2.15) and (2.5). The above transition H{a) — > H{a) generates a new operator H{a) as 
compared to H{a), but Eqs. (2.17) and (2.18) show that W{a)H{a)W^{a) turns out to 
be equal to the original total Hamiltonian (cf. [p5||). 

We would like to stress that the transformation WHW^ under consideration should 
not be understood here as W{a)H{a)W^{a). The latter would be a new operator H'{a), 
which, in general, does not coincide with H. For a detailed discussion of different unitary 
transformations, see Sect. 6. 

"^In cumbersome formulae summations over the dummy spin indices are sometimes omitted 
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2-4 Elimination of "Bad" Terms 

The next step is to fulfil the requirements i) - ii). If we want the no-clothed-particle 
state fl and clothed one-particle states to be H eigenvectors, the r.h.s. of Eq. (2.19) must 
not contain some undesirable terms. Particularly, K{a) must not contain the b^^dc^ac^- 
type terms because they would give rise to the ffm state^, when acting on Q, and K{a)Q 
could not be proportional to Q. Similarly, bjbcac^ converts a one-fermion state bjQ into 
a fm state. But just terms of this kind enter into the operator V{a) which occurs in the 
r.h.s. of Eq. (2.19). In this connection, recall that Hj{a) = V{a) + Mren(tt) where V{a) 
is derived from V{a), see Eq. (2.16), by means of the replacement a ^ a. 

As we have argued above, the terms 

b^b^aj, bMc^cic, bMcdc, dcdjaj (2.20) 

in V{a) do not allow the clothed no-particle and one-particle states to be H eigenvectors. 
The remaining terms in V{a) are Hermitian conjugate of (2.20). We shall call "bad" all 
these terms. The contribution bc^dc^ac^ will be called the "bad" term of the class [3.0]: it is 
a product of three creation operators with destruction operators not included. The other 
three terms in (2.20) belong to the class [2.1]: two creation operators and one destruction 
operator. 

The interaction Hj includes also the mass counterterm M^en (see Eq.(2.7)). The latter 
contains bad terms of the class [2.0] (see,e.g., Eq.(2.27)). The self-energy correction to 
the particle mass can be represented by a series which starts with the terms of the 
-order(see Eq. (2.31)). So, Mren{a) ~ 9^ while V{a) ~ g^. 

Let us eliminate from K[a) the bad terms of the g^ -order. For this purpose we choose 
such R that 

V + [R, Hf]=0 . (2.21) 

One readily verifies that Eq. (2.21) cannot be satisfied until R{a) is linear or bilinear 
in a. To meet the equation, R{a) must be a three-operator, e.g., have the structure of 
V{a). Then, the commutator [R,Hp] will also be three-operator expression since Hpia) 
is the two-operator 0. 

Let us assume that the antihermitian R{a) contains the bad terms of the same kind 
as V{a) is. Namely, we put R{a) = 1Z — V) where (cf. Eq. (2.16)) 

n = j dp'dpdk Y.{R^^{p'r';pr)bJ{p',r')b,{p,r) + R^^{pV;pr)b,^{p',r')dJ{-p,r) + 

r'r 

^with a transparent abbreviation ffm for a "fcrmion-antifermion-meson" . 

^Note that the commutator of the m-operator term and n-operator one yields a (m + n — 2)-operator 
contribution 
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+ R^-,{p'r';pr)d,{-p',r')b,{p,r) + Rg{p'r';pr)dc{-p\r')dJ{-p,r)} a,{k) . (2.22) 

The c-number coefficients Rfj{i,j — 1, 2) are to be derived from Eq. (2.21), see Appendix 
A. We find tfiat the solution exists if /i < 2m. This condition has a clear physical 
meaning, viz., the meson can decay into the ff-pair if // > 2m, and, therefore, one-meson 
state cannot be stable,i.e., it cannot be an H eigenvector. Once [R, Hp] — —V, Eq. (2.19) 
can be rewritten as 

K{a) = Hpia) + M,e„(a) + ^[R, V] + [R, Mren] + ^[i?, [R, V]] + ... . (2.23) 

Thus we have removed from K{a) all the bad terms of the g^-order. 

However, the r.h.s. of Eq. (2.23) embodies other bad terms of the g"^- and higher 
orders. For example, [-R, ^] contains the terms g'^bc^dc^ac^aj of the class [4.0], which do 
not destroy the physical vacuum fl (our evaluation of [R, V] is given in Appendix A). In 
addition, we find in [R, V] the terms g'^bc^dc^ac^ac of the class [3.1], which neither destroy 
ajQ, nor retain it with a multiplicative factor. These and similar bad terms can be 
eliminated in a way analogous to the described above via one more transformation 

ap = W4{a')apWj{a') , (2.24) 

where W4 = cxpi?4 and R4 is an expression of the 5'^-ordcr, which consists of the above 
bad four-operator terms. This R4 should be such that [R4, Hp] would cancel the latter 
terms. 

Note also bad terms of the classes [2.0] and [1.1], which are present in M (see Eq. 
(2.27)), and similar terms, which appear after normal ordering of [-R, V^] (see Appendix 
A). The bad terms of the class [2.0 ] must be removed from K[a) as well. We shall show 
in Subsect. 2.5 how they may be cancelled under a condition that relates the physical 
masses with the input parameters mo, fJ>o and g. 

Further, the double commutator [R, [R, V]] in Eq. (2.23) is composed of five-operator 
terms (cf. the footnote 3), and there are bad terms among them. In particular, after 
reshuffling the operators into normal order, new three-operator bad terms occur. However, 
they are of the g'^-order. This type of bad terms can also be found in [R, M] . The 
subsequent unitary transformation makes it possible to remove all bad terms of the g^- 
order. 

Along the guideline, one may eliminate from the Hamiltonian the bad terms of increas- 
ing orders in the coupling constant g. It is assumed that in the limit the requirements i) 
and ii ), which are equivalent to the absence of bad terms in K, will be fulfilled. 

Finally, if our clothing procedure were perfect, the resulting representation K of the 
total Hamiltonian would possess the property 

K{a)\k)^ = Hp{a)\k)^ = Uk\k)^ (2.25) 
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with I k)^ = aJ{k)Q . In other words, the new interaction term Ki{a) = K{a) — Hpia) 
would satisfy the equation Kj{a) \ k)^ = 0. Analogous equations will hold for the physical 
vacuum Q and the clothed one-fermion and one-antifermion states. 

2.5 Particle Mass Renormalization 

The cancellation of the bad two-operator terms in K{a) will be demonstrated for 
those of them which are bilinear in meson operators Oc and a^^ . The latter originate , first 
of all, from the meson mass counterterm 

M^es = - Ai') / <p\x)dx (2.26) 

Indeed, substituting the expansion (2.9) for (j){x) into Eq.(2.26), we obtain 

2 2 

Mrnes = I ^° ~^ [2a,^{k)a,{k) + a,{k)a,{-k) + a,t(fc)a,t(-fc)]rffc (2.27) 

As shown in Appendix A, terms of the same operator structure occur in K{a) after normal 
ordering of the commutator [i?, V] from Eq.(2.23): 

H 
'k 



f ^[2aJ {k)ac{k) + ac{k)ac{-k) + aj {k)aj {-k)]dk , (2.28) 



where is determined by Eq.(A.20). 

Now, we see that the sum of (2.27) and (2.28) gets equal to zero if 

/io^ - = -4tfc or fi^ = /io^ + (2.29) 

Here the quantity tf^ must not depend on k along with /i^ (cf. our argumentation in 
Subsect. 2.2). Unfortunately, this independence is not automatically provided in the 
Schroedinger picture that we use throughout. The integral that determines tf^ (see 
Eq.(A.20)) is quadratically divergent and one needs to overcome the trouble (e.g., by 
introducing a cutoff factor). So, special efforts are required to yield the proof of indepen- 
dence (see, e.g., Ch.6). 

In the same way, the fermion mass countrerterm 



Mferm = ("^0 — rn) J 'ip{x)^jj{x)dx (2.30) 

cancels, under a proper condition, all the terms bilinear in the fermion operators, which 
arise from | [R, V] result of normal ordering. 

Up to now we have considered the bilinear terms of the (yf^-order. Normal ordering of 
the six-operator and other terms of K{a) gives bad two-operator terms of the g"^- and 
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higher orders. To ehminate them we suppose that = /iq^ — /i^ and 6m = mo — m may 
be expanded in the series 

oo oo 

^/i' = E ^?'"('^/^')2n, Sm=Y^ ^7'"(5m)2„ (2.31) 

n=l n=l 

Let us assume that (5/i^)2 and (5m)2 are used to remove the two-operator terms of the 
(jf^-order as described above. Then, the terms (5yU^)4 and {Sm)^ are destined to cancel 
two-operator terms of the ^f^-order, and so on. 

2. 6 Some Remarks 

So, the transformation reahzed hy W = expR fulfills an incomplete clothing, viz., it 
removes bad terms of the least order in g. The no-particle fl and one-particle states bjQ, 
djQ and ac^Q constructed at this stage are merely approximate H eigenvectors. 

We shall confine ourselves to the consideration of this transformation only while dis- 
cussing the bound state problem like the deuteron (Sect. 4). Even this simplest applica- 
tion of our approach turns out to be rather cumbersome. 

This section has been aimed to show how the unitary transformations of the original 
Hamiltonian for a system of interacting fields can be regarded as the introduction of new 
creation(destruction) operators instead of the initial "bare" ones. These new operators 
and corresponding one-particle states occurring at the first stage of the clothing procedure 
may be called "partially clothed". We juxtapose the one-particle states to the observable 
particles (say, pion and nucleon). 

Let us note that the clothed states (operators) are not the in(out)-states (operators) 
of RQFT (see, e.g., [|^, Ch. 17). Indeed, the two-particle in-states ain"^ {ki)aiJ {k2)fl are 



H eigenstates while the two-particle clothed states aJ{ki)aJ{k2)Q do not (simple models 
with noninteracting particles are the evident exception). Within the in(out)-formalism 
it is supposed that H has the following eigenstates, viz., no-particle state (physical vac- 
uum), states with one in-particle, two-, three-, etc., these states being analogous to the 
corresponding Hq eigenstates. Meanwhile, the clothing formalism needs not such a sup- 
position: one can find explicit expression for no- and one-particle clothed states in terms 
of bare states (using formulae of this section under the condition /i < 2m ). In addition, 
the in(out)-formalism does not consider two-particle states which have no Hq bare part- 
ners, e.g., an in-state describing the deuteron. In the clothing formalism the problem of 
deuteron-like states is subject to further investigation (see Sect. 4). 

3. Generators for Space Translations and Space-Time Rotations 
within the Clothing Procedure 
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In the previous section we have expressed the total Hamiltonian H in terms of the 
clothed operators. H is the time translation generator of the Poincare group. Here we 
shall discuss how the rest of the group generators (total linear and angular momenta, 
and generators of Lorentz boosts) depend upon the clothed operators. This will allow 
us to formulate the transformation properties of the clothed no-particle and one-particle 
states under the Lorentz boosts. First of all, let us consider various constraints imposed 
on clothing transformations by the general symmetries in RQFT and those which are 
specific for a given field model. 

3. 1 Total Momenta and Other Motion Integrals in Terms of Clothed Operators 

Similar to the determination of the total Hamiltonian if as a function K{a) of the 
clothed operators a we can obtain the expression Pc{oi) for the total linear momentum 
P (cf. Eq.(2.17)): 

P(a) = P{W{a)aW\a)) = W{a)P{a)W^{a) = Pc(a) (3.1) 

One can show that WPW^ is simply equal to P. Indeed, this is equivalent to WP = PW 
that follows , in its turn, from [R, P] = (remind that W = expR). The validity of the 
latter can be verified using Eq.(2.22), (A.l), (A. 2) and the well-known expression of the 
total linear momentum in terms of the creation(destruction) operators 

P = J ka^{k)a{k)dk + j p^P {P,r)h{p,r) + d\p,r)d{p,r)\dp (3.2) 

(see, e.g., Eq.(7.33) in ||15|). Alternatively, one can take the representation (A. 4) for R 
and use the equation [P, V^] = that holds because of the invariance of V with respect 
to space translations. 

So, we have found that Pc(«) = P(a). This means that P is the same function of the 
clothed operators as of the bare ones. 

Analogous statements are valid for the total angular momentum M [] and the 
baryon(fermion) number operator B = J il)\x)il){x)dx. This means that the clothed 
states f2, Oc^fi, hc'^VL, and dc'^VL have the following properties: a) they are eigenvectors of 
P ; b) they are transformed under space rotations in the same manner as the relevant 
bare states do ; c) they possess definite B values. 

Instead of verifying the properties [R, P] = [R, M] = [R, B] = with the solution 
R of the equation [R, Hp] + V = one may consider them as some new requirements 
supplementary to those listed in Subsect. 2.1 as i) - iv). These requirements would 

^The field-theoretical formula for the generator can be found in [^(Ch.ll), see also Subsect. 3. 2. 
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result in definite restrictions on the coefficients involved in the expansion (2.22) for R. 
For example, it follows from [-R, -P] = that R^{p'r';pr) must have the form 5{p + 
k — p')r^{p'r';pr). The equation [R, M] = means that Rf^{p'r';pr) must depend on 
rotationally invariant combinations of its arguments. Besides the condition [R, B] = 
prevents R to be dependent on terms of the hc^dc- kind (we eliminate such terms from 
the beginning assuming the form (2.22) for R). 

One may add to these restrictions those which are consequences of the following re- 
quirements: clothed operators and clothed one-particle states must have the same trans- 
formation properties with respect to space inversion, time reversal and charge conjugation 
as their bare partners. 

Let us stress that all the constraints are exact whereas the equation [i?. Hp] + = 
considered in Subsect. 2.4 is merely approximate one. However, its solution has all the 
properties in question since the interaction V commutes with P, iW, B, etc. (see Eq.(A.4)). 

3.2 Transformations of Bare Operators and States under Lorentz Boosts 

A distinctive feature of the conventional relativistic dynamics ("instant" form af- 
ter Dirac ) is that the generators N = {N^,N'^,N^) of the Lorentz boost A = 
exp{i(3N)^ contain interaction terms while the linear P = (P^, P^, P^) and angular 
M = (M^, M^, M^) momenta are determined by the same expressions as for free fields. 

In order to see this explicitly let us resort to the Lagrangian formalism where the quan- 
tities of fundamental importance are the energy- momentum density tensor T^'^(x) and 
the angular momentum density tensor 7W'^'^'^(x)| (see, e.g., Ch.ll in [Q). By definition, 
= {H,P) = Jr^^'{x)dx, = ejkiM''\ and = M°^' where M^"' = J M'^f'^{x)dx. 
According to the Noether theorem all the ten operators are time independent, i.e., they 
are the motion integrals. In other words, they can be evaluated at t = 0, i.e., they can be 
expressed through field operators in the Schrodinger picture. 

The corresponding representation of N depends on the form of T^'^(x) (nonsym- 
metrized or symmetrized]^, that is utilized. Here we shall use the nonsymmetrized form 
(see Eqs. (13.45) and (13.47) in ||38|), which leads to 




^Here (3 ~ l3n, n — ^ and thp = v, where v is the velocity of a reference frame moving along the n 
direction. In this paper we use the system of units in which the light velocity c is equal to unity. 
^Greek labels run the values 0,1,2,3. 

^°The symmetrized form with the Belinfante ansatz for Lorentz boosts has been employed in our 
talks ^ . Another application of the form can be found in a covariant description of electromagnetic 
interactions with nuclei pO| . 





(3.3) 
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where Np is the free part of N: 

NF^Nf,rm + N^^, , (3.4) 

Nferm — — j X'tl!{x)[—ij'V + 'm\ijj{x)dx + ^ j '4^{x)'^'4^{x)dx (3.5) 

Nmes ^ [^'(^) + ^(t>{xf + /uV'(a;)] dx (3.6) 

In accordance with the particle mass renormahzation described in Subsect.2.5 we have 

separated the contribution to N from meson and fermion mass counterterms (cf. Eqs. 
(2.26) and (2.30)): 

iV,,„ = AT— + Nilir (3.7) 

NT:: = -\W - / xct>\x)dx (3.8) 

NlZ;^ = -(mo - m) y xij{x)ilj{x)dx (3.9) 

Now, to express the generators through the creation(destruction) operators let us take 
the expansions (2.9) - (2.11) and employ the relation 

j x^ex]){iqx)dx = — i(27r) -Q—S{q) 

Then, e.g., we find 

NL.^y.HyM,fJ^^^^^Pi,'-,W,, (3.10) 



or 



.da\k) ,,,,5a(/e). 



NLes - \ j -ki^-ik) - aKkf-^]dk (3.1O0 
Simultaneously, we get (cf. Eq. (2.16)) 

Nj = -J x^V(x)dx = - 

Q f Tfl d 

(2^ J {2ukEp,Epy/- dkJ^^P + ^ - ^'^ Mp'r'hMpr) b\p'r')b{pr) + 

+ u{p'r')'^^v{—pr) h\p'r')d\—pr) + v{—p'r')^^u{pr) d{—p'r')b{pr) + 

- v(-p'r')-f5v(-pr) d^(-pr)d(-p'r')} [a(k) + a^(-k)] (3.11) 

with the Yukawa interaction density V{x) (see Eq. (2.2)). 

These formulae enable to perform directly transformations of the bare operators and 
states under the Lorentz boosts. In particular, in the infinitisemal case with | |<C 1 (j 
= 1,2,3) one has 

AOo = exp[i/3{NF + Ni + Nren)Po 
16 



~ [1 + lf3{Ni + Nren)]^^ ^ (1 + if^N 1)^0 , (3.12) 

where we took advantage of the relation 

Nf^o = , (3.13) 

and omit terms of the order higher than g^. 

Eq. (3.12) means that the bare vacuum is not invariant with respect to A, viz., 
Af^o 7^ ^0- A moving observer "sees" AQq as the superposition of no - particle state Qq, 
NNn states, etc. 

Similarly, bare one-particle states (e.g.,a^(A;)no) are not transformed with respect to 
A as in the free case where along with Eq. (3.13) one has the property Q 

e'^^^a^{k)e''f^^^ = a\Lk) (3.14) 

In Eq. (3.14) L denotes a pure Lorentz transformation ("boost") with the matrix: 



M° : Uj 



(3.15) 



3 l+uo 

where u'^ = (m°, u) = {chjS, nsh(3) is the four - velocity vector. The boost converts the 
four - momentum k = {uJf^, k) into k' = Lk = (t^^', k'). 

We have Aa'^{k)Qo = aA^(A;)Af2o with the transformed meson operator aA"''(fc) = 
Aa^(fc)A"''. For an infinitisemal boost the operator 

a^^{k) = e'^^ a\k)e-'^^ ~ a^k) + il3[Nmes. a\k)] + i(3[Ni, a\k)] (3.16) 

contains bare fermion operators due to the last ("interaction") term in the r.h.s. of 
the equation. By using Eq. (3.10') one can see that the two first terms yield a^{k) — 

f^^^^dff^^k — '^^(^k ~ ^ ^ '^fc'^) ^^^^ coincides with the operator a\Lk) in Eq.(3.14) 
for the infinitisemal L. 

Now, taking into account Eqs. (3.12) and (3.16), one can ascertain that the trans- 
formed state Ka\k)VLQ is composed of the one - meson state a^{k')VLQ with the properly 
changed momentum k' = k~ ujj^/S and the states | //) and | //tttt) containing fermions. 

3.3 Boost Generators for Clothed Particles. Elimination of Bad Terms 



^^Under the discussion it is convenient to proceed with the operators a(k) which obey the covariant 
commutation relation [a(fc), a^(fc')] = ujj^S{k — k') (cf., e.g., Eq. (7.23) in [|l^). In the framework of our 
consideration it is equivalent to replacement of a{k) by a{k) / ^JUJ^. It is true for the respective clothed 
operators, so that, for instance, the first of the relations (2.5) should be replaced by [ac(fc), ac^(fc')] = 
uJi^Sik-k'). 
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It is reasonable to anticipate that the physical vacuum and clothed one-particle 
states (e.g., aJ{k)Q) should be, respectively, the no-clothed-particle state and clothed 

one-particle states from the point of view of a moving observer. More exactly, they 
should meet the relations 

An = n (3.17) 

and 

AaJik)n^aJ{Lk)n (3.18) 

The previous experience of handling with Eqs. (3.12) and (3.16) prompts that these 
conditions could be provided (at least, approximately) if we shall manage to remove 
bad terms from N (in practice, some of them) while expressing it through the clothed 
operators. In this connection, let us write down 

N = N{a) = WN{a)W^ = B{a) = e^^^^Npia) + Ni{a) + Ar^en(a)]e"^^"^ 

^Nf + Ni + [R,Nf] + [R,Ni] + ... (3.19) 

(of. Eq. (2.19)) and then remove from the r.h.s. of (3.19) all the bad terms of the 51^-order 
by requiring that 

[Nj,, R]^n{^- j x^V{x)dx (j = 1, 2, 3) (3.20) 

Now, we want to show that Eqs. (3.20) will automatically hold if R satisfies the 
condition (2.21). For this purpose, one can use the representation (A. 4) for such R: 



dte-'^ J V{x)dx , (A4) 



R — —i lim 

e->0+ Jo 

where V{x) — V{x,t) — e^^^^V{x)e~^^''* is the interaction operator in the Dirac picture. 
Being a scalar it is transformed as 

under the Lorentz boost L determined by the matrix (3.15). Note that one can directly 
verify the validity of Eq. (3.21) in case when Np is the Schroedinger operator (as it does 
in Eq. (3.19)) while V{x) being any Lorentz invariant operator in the Dirac picture. 
In order to exploit this property note that 

[N'f, R] - -^^K^^i?e-''^^^] 1,,^, (3.22) 

or taking into account firstly Eq.(A.4) and then Eq.(3.21), 

[Nl,R] = - lim lim — / dte~'^ / V{Lx)dx (3.23) 
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For the infinitesimal boost we have 

V(Lx) = V(x -f3t,t- (3x) = V(x^ - pH, x^ x^ t - (3^x^) (3.24) 

and 

[ATI , R\^- hm / dte-'^ / hm ^Vix^ - (3H, a;^ x^ i - /?^xMda; 
e->o+ Jo J /?i^o ap 

//•oo f) /•oo /• 

a;Ma; / e~'^—V(x,t)dt+ tdte''^ —Vix,t)dx\ (3.25) 
io at Jo J oxi 

One can show (cf., the proof of the relation (A. 4)) that the first term in the square 
brackets yields — / x^V{x)dx. At the same time the second term is equal to zero since 
the operator V{x,t) (more exactly, its matrix elements) vanishes at x^ — ±oo. So, we 
get the desirable relation (3.20) with j — 1 (the cases j = 2, 3 are analogous). 

So, the transformation W — expR eliminates simultaneously the three-operator terms 
~ both from the total Hamiltonian K{a) and the boost generators B{a). One should 
emphasize that the proof is valid for any Lorentz scalar function V{x, t). Specific expres- 
sions for Np have not been required as well since all we have needed is Eq. (3.21). 

After this elimination of bad terms we get by analogy to Eq. (2.23), 

B{a) = NF{a) + N^enia) + ^[R, Ni] + [R, Nren] + l[R, [R, Ni]] + ... (3.26) 

We shall not exemplify separate interaction terms in the r.h.s. of this equation since their 
structure repeats that for the corresponding contributions to K{a) (cf., e.g., Eqs.(3.11) 
and (2.16)). 

4. Equations for Bound and Scattering States in RQFT 

The clothed one-particle states are eigenstates of the the total Hamiltonian H ac- 
cording to their definitions (see Subsect. 2.1). There may be other H eigenstates which 
describe physical systems resembling one-particle states, viz., the states with discrete 
values of the system mass that may be defined as the system energy in the rest frame of 

reference (scattering states belong to continuous values of the mass). 

First of all, we keep in mind the simplest bound states similar to the hydrogen atom or 
the deuteron. In the non-relativistic approach the wavefunction of such two-body state 
is the product of a function that describes the system as a whole and other function being 
dependent on the internal variables (for instance, the relative momentun p = j{Pi — P2) 
for two identical particles). Therefore, the ccntre-of-mass motion is separated from the 
internal motion. It is not the case for any relativistic model which satisfies the Poincare 
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algebra . In fact, the coupling between the internal and centre-of-mass motions is inherent 



to relativistic theories of interacting particles (see, e.g., the papers [^, where within 
simple field models the two-body bound states are studied in a moving reference frame). 

For the Yukawa model the corresponding states may be fermion-fermion states 
(deuteron-like), meson-fermion ones, etc. Of course, we are not able to find the ex- 
act Hamiltonian eigenstates, except some exactly solvable models (see Q)- However, 
within the clothing procedure in question one can suggest reasonable approximations to 
this problem. 

4.1 New Zeroth Approximation for the Total Hamiltonian K 

Our approach is based on the choice of an appropriate zeroth approximation (ZA) to 
the total Hamiltonian expressed through the clothed operators, i.e., the operator K de- 
termined by Eq. (2.17). Since its free part K2 = Hp{a) has no deuteron-like eigenstates, 
we shall try to take 

KzA = K2 + g^Kf^ (4.1) 
by adding to the two-operator(one-body) contribution 

^2 = y" uj]^aj{k)ac{k)dk + J EpJ2[bJiP, r)bcip, r) + dj{p, r)dc{p, r)]dp = + 

(4.2) 

the four-operator (two-body) contributions of the - order which arise from the com- 
mutator |[_R, V] = ^[Rs, V] in the r.h.s. of Eq. (2.23). This commutator is evaluated in 
Appendix A. Doing so, we obtain the decomposition 

= g^K^ = K{NN NN) + K{NN NN) + K{NN NN) 

+K{ttN ttN) + K{ttN ttN) + K{7rn NN) + K{NN rnr) (4.3) 

with the separate interactions between the different clothed particles. They are displayed 
(very schematically) in Fig.l where the graph (a) represents the nucleon-nucleon interac- 
tion 

K{NN NN) = J2j dp[dp2dp^dp2 

r,r' 

Vnn{p'i, r[,p'2, r2;Pi,ri,p2, r2)bj{p[,r[)bj{p'2, r2)6c(Pi, ri)6c(P2> ^2) (4.4) 
while the pion-nucleon interaction 

K{ttN ^nN) = Y,J dk'dp'dkdp 

r,r' 

V^k', p', r'; k, p, r)aj{k')bj{p', r')a,{k)b,{p, r) (4.5) 
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is displayed by the graph (c). 

Exphcit expressions for the coefficients Vt^n and Vnn will be given in the next Subsects. 

(2) 

Here, however, one should note that all these terms of K^^ describe only real processes 
such asiV + iV^A^ + A^, 7r + A^->7r + A^, etc. The bad terms of [R3, V] are not 
included in Kza for the reasons discussed in Subsects. 2.4 and 2.5. For example, the 
terms of the kind [4.0] (e.g.,b'^d^b'^d^ and b'^Sa^a'^) and [3.1] (e.g., b'^Sb'^b and b'^d)a)a) 
must be removed during our clothing procedure via the transformation W4 = expi?4 (see 
Eq. (2.24)). Recall that [i?4,i^'2] niust cancel these bad terms. 
In its turn, the transformation W4 (see Eq. (2.23) for K) 

W^K{a')W^^ = W4{K2{a') + ^[R^ia'), V{a')] + ...)1^4^ = 

K^ia') + ^[Rsia'), V{a% + [R,, K2] + ^[R^, [R3, V],] + ... , (4.6) 

brings in the Hamiltonian new four-operator terms in addition to those mentioned above, 
see Eq. (4.3). Here, [R3, V]^ denotes the four-operator part of [R3, V]^ After the normal 
ordering, the double commutator [R^, [R^, V]^] yields new four-operator interactions in 
the total Hamiltonian. However, they are of the g"^ - order whereas the interaction terms 
which we have included in Kza are of the - order. So, the latter are not altered by W4. 
They exhaust all the interaction terms of the g^ - order, which remain in Kza- 

We hope that Kza eigenstates are good approximations to exact K eigenstates. The 
former should be found nonperturbatively (e.g., by means of numerical methods). Five- 
operator and more complicated interaction terms can be taken into account via pertur- 
bation theory recipes. 

4.2 Meson-Nudeon Eigenstates of Kza cind Pion-Nucleon Quasipotential 

The operator Kza has an important property: it conserves the total number of clothed 
particles. In particular, Kza transforms clothed two-particle states (e.g., of the A^A^ or 
TiN types) to two-particle ones. Moreover, the Fock subspace of all the clothed states can 
be divided into several sectors (the A^A^ sector, the vrA^ sector, etc.) such that Kza leaves 
each of them to be invariant, i.e., for any state vector $ of such sector Kza^ belongs to 
the same sector. 

Let us show in the simplest case of the nN sector that the property just mentioned 

allows us to reduce the eigenstate equation Kza^^ = to the related Schroedinger 

^^The two-operator terms of [i?3,y] are supposed to be cancelled with the respective mass counter 
terms. 
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equation of the particle-number-conserving quantum mechanics. For this purpose we 
seek <l>^ as the following superposition of the nN sector states Oc^&c^^, 

^nN = J2 f dkdp^^j,{k;p,r)aj{k)bj{p,r)n (4.7) 

In the vrA^ sector Kza is equal to + + K{ttN vrA^) because the rest terms of 
KzA give zero when acting on Thus, the equation Kza^^n — ^^wN reduces to 

[K^ + Kr, + K{7fN 7riV)]<l>f^ = (4.8) 

Taking the scalar products of the both parts of (4.8) with {aj {k)bj {p,r)Q |, we get 
the relevant equation for $^^(fc;p, r), 

{E-u;}^-Ep)<l>^^{k;p,r)= J] J dk'dp'VMk,P,r; k' ,p' ,r')<i>^j,{k';p' ,r') (4.9) 

r' 

The kernel V-j^n of this integral equation is determined by Eq. (4.5) (an explicit expres- 
sion for it is given below). This kernel can be called a quasipotential: in the coordinate 
representation it may depend not only on the particle coordinates but on their derivatives 
as well. In our opinion, a popular name "effective Hamiltonian" is inappropriate for Kza- 
When one deals with an effective Hamiltonian one has to argue that its eigenvalues and 
eigenvectors coincide (at least, approximately) with those of an original Hamiltonian. In 
the framework of our approach we do not need such a proof since we believe that Kza is 
in a sense a major part of the total Hamiltonian H and, therefore, the above approximate 
coincidence is provided. 

The solutions of Eq. (4.9), which belong to a discrete spectrum (if it exists), describe 
nN bound states. As emphasized above, the solutions should be found nonperturba- 
tively. Continuous vrA^-mass values correspond to vrA^-scattering states. The respective 
S'-matrix elements may be evaluated either exactly by using numerical methods of solving 
the Lippmann-Schwinger equation for the T-matrix with the interaction K{txN — » vrA^) 
(see, e.g., [0), or approximately in the framework of old-fashioned noncovariant pertur- 



bation theory (see, e.g., Ch.l in |^ or Ch.ll in 

To obtain the explicit expression for K{itN ttN) , one needs to seperate out the 
ajbjacbc-kind terms of the commutators [7?., V"!^] and [7^, V"'^]^ (see Eq. (A. 14)) 

K(7iN vrAT) = i : |[7^, V^].iv + [7^, V^ll^j : 

= — ^ y dk2dp2dkidpi 

!^[V-^\R%^{p2r2;p,r,) + [V-^\R^%{p2r2;p,r,)yj{k^^^ , 

(4.10) 
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where the symbol : : denotes the normal orderin^^. 

This result has been obtained in |T^(cf. Eq.(A.7) therein). The corresponding co- 
efficients VirNik' ,p' ,r'] k,p,r) that determine the pion-nucleon quasipotential (see Eq. 
(4.5)) are equal to 

{aj{k')bj{p', r')n I K{nN ^ vriV) | a,{k)h,{p, r)fi) (4.11) 
They can be represented in the following covariant (Feynman-like) form 

1 Til/ 

VTrN{k2,P2,r2] ri) = 777^^(fc2 + P2 - ^1 -Pi) 



ri r 1 

u{P2,r2)i 



1 

+ 



}^2+^2 + 'm }^i+^i+m 



1 

+ 2 



1 

+ 



2-//i + m ^1-^2 + ^ 



(4.12) 

(see also Appendix in P^). 

In order to comment this expression, let us consider the Feynman graphs in Fig. 2 for 
the S'-matrix elements of vrA^ scattering. According to Feynman rules the four-momentum 
of the internal nucleon line in the graph 2a equals either the sum ki + pi of the incoming 
four-momenta or the sum k2 + P2 of the outgoing four-momenta. These sums are equal 
due to the energy and momentum conservation, viz., the S'-matrix contains the multiplier 

5{k2+P2 - ki -pi) = 5{ujf^^ + Ep^ -^ki ~ Ep^)^{k2+P2 - fci - Pi) (4.13) 

Therefore, the Feynman propagator corresponding to the internal line can be written 
either as (Jfi -t-^i + m)~^ or as (^2 + ^^2 + ^«)~^. In the case of quasipotential the energy 
conservation is not assumed (only the total three-momentun is conserved) and hence 
ki + pi is not necessarily equal to ^2 + P2- The representation (4.12) shows that V^tv 
includes the contribution associated with graph 2a and it can be obtained if we juxtapose 



to the internal nucleon line the half-sum | 



+ 



In the case of S'-matrix we juxtapose to the internal line in graph 2b either the 
propagator (^2— ^i+"^)^^ or (^i— ^2+"^)^^. In the case of quasipotential pi—/c2 7^ P2~ki, 
in general, and the half-sum of these propagators does correspond to the internal line (it 
is the second half-sum in the curly brackets in the r.h.s. of Eq. (4.12)). 

It follows from these observations that multiplying Ktat by the factor —2'ki5{uj^^ + 
Ep^ — ojj^^ — Ep^ ) we shall obtain the S'-matrix elements for TriV scattering in the g"^- 
order. 

So, we have seen that the r.h.s. of Eq. (4.12) resembles the Feynman amplitudes being 
different from them in the two respects: i) the multiplier 5{k2+P2 — ki—pi) is substituted 
^•^Henceforth summation over the dummy spin indices is implied. 
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instead of (4.13); ii) the above Feynman propagators are replaced by the corresponding 
half-sums. 



4-3 Clothed Nucleon-Nucleon Eigenstates and Nucleon-Nucleon Quasipotential 

Now, we consider the Kza eigenstates which belong to the A^A^ sector, being super- 
positions of the kind, 

^NN = J2 dp-^dp^^NNiPi, ri; Pa, r2)bJiPi,n)bJ{p2, rs)!! (4.14) 

A subset of such states with a definite momentum P is determined by Eq. (4.14) with 
the coefficients <l'ArAr(l;2) ~ 5{pi + P2 — P) \^ Obviously, Kza^nn is the the state 
vector of the same sector. In fact, the operators K^^, K{ttN — * vrA^), K{ttN vriV), 
K{NN NN), K{7nT NN) and K{NN tttt) involved in Kza do not contribute 
to KzA^NN (see Eqs. (4.1) and (4.3)), and we find that Kza is reduced to the operator 
Kn + K{NN — > A^A^). So, the eigenvalue equation Kza^nn = ^^nn yields the equation 



Kn + K{NN ^ A^A^) 



= (4.15) 



in the sector. 

The corresponding equation for ^ffj^{l;2) (see Eq. (4.14)) can be derived 
from Eq.(4.15) by taking scalar products of the both parts of the latter with 
{bJ{Pi,ri)bJ{p2,r2)Q |. Doing so, we get 

{E - Ep^ - EpJ<l>f ^(pi, n; r2) 

= J dp[dp'2VNNiPi,ri,P2,r2;p[,r[,p'2,r'2)(^f;^{p[,r[;p'2,r'2) (4.16) 
with the properly symmetrized interaction (the quasipotential) 

Vnn{1, 2; 1', 2') = ^V^Nil, 2; 1', 2')-V^7v(l, 2; 2', l')-Viv^(2, 1; 1', 2')+Viv^(2, 1; 2', 1') 

(4.17) 

for the two clothed nucleons. 

The two-body operator K{NN — > A^A^) (see Eq. (4.4)) is generated by the second 
term in the curly brackets of Eq. (A. 14) and its H.c. : 

K{NN ^NN) = ^: |[7^, VV^v + [n,V%^^ : (4.18) 

One should note that the coefficients VAr7v(l, 2; 1', 2') in Eq. (4.4) are not in the one- 
to-one correspondence with K{NN —>■ NN), viz., they can be changed without altering 
the latter. For instance, the property 

6(1)6(2) = -6(2)6(1) (4.19) 



*Here and sometimes below we use the evident abbreviations, viz., 1 = {pi,ri), etc. 
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enables one to replace VArAr(l, 2; 1', 2') in Eq. (4.4) by —Vnn{^, 2; 2', 1'), and so on. 

Moreover, the operator K{NN NN) remains unaltered when adding to Vat at 
arbitrary functions ^^(l, 2; 1', 2') or Sfi{l,2;l' ,2') which are symmetrical under the 
transpositions 1 ^ 2 or 1' 2', respectively. The above mentioned replacement 
VArAr(l, 2; 1', 2') — Vat at (1, 2; 2', 1') is equivalent to the replacement 

Vnn{1, 2; 1', 2') ^ Viv^(l, 2; 1', 2') + Sn{l, 2; 1', 2') (4.20) 

with Sail, 2; 1', 2') = -Vnn{1, 2; 1', 2') - Viv^(l, 2; 2', 1'). 

A distinctive feature of the coefficient (4.17) is its invariance with respect to the 
transformation (4.20) with arbitrary Sji. 

After these notations, we write down one of the possible expressions for Vnn, that can 
be obtained using Eq. (4.18), 



Vnn{1,2;1',2') = 

1 I ^ ^ + F ^ -]vnHp.rup[r[)v'f 

2 J [Ep^- Ep>^ - uj^ Ep>^ - hp^ - uj^ J 

The respective quasipotential is 

y^^(l, 2; 1', 2') = --l—Mv^ +P2- p[ - P'2) 



(27r)3 2JEp^Ep^Ep^,Ep^; 



Pi P2 Pi' P 

-(1 ^ 2) (4.22) 

Expression (4.22) is the A^A^ part of an one-boson-exchange interaction derived via the 
Okubo transformation method in |^ (cf. [0). The potential Vtvat consists of the direct 
term written explicitly and the exchange term (1^2). In order to obtain the latter one 
needs to replace Pi, ri by P2, r2 and P2, r2 by p^, ri in the former. 

As has been pointed out in [005 a distinctive feature of the potential is appearance 
of a covariant (Feynman-like) "propagator" 

1 r 1 1 1 

+ ... ... > (4-23) 



2 [ (p'l - pi)2 - /i2 (p'2 - P2Y - 

where p = [Ep, p) is the nucleon four-momentum. On the energy shell that is when 

E, = Ep^ + Ep^ = Ep,^ + Ep,^ = Ef (4.24) 

^^There one can find another representation of the nucleon-nucleon quasipotential, which resembles 
the expressions of old-fashioned perturbation theory (see, e.g., Ch. 13 in |p5[). 
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the r.h.s. of Eq. (4.23) becomes the genuine Feynman propagator which appears when 
evaluating the S-matrix for A^A^ scattering in the - order. The respective graphs 
are displayed in Fig. 3. Like Kta? the quasipotential Vnn can be associated with these 
Feynman graphs being different from the corresponding Feynman amplitude in the two 
respects, viz., Vnn does not contain 5{Ep^ + Ep^ — Ep^ — Ep'^), and "propagator" (4.23) 
now corresponds to the internal meson line in graph 3a. A more extended analysis of this 
observation has been been given in [p5|] . 



In conclusion, one should note that all the quasipotentials are nonlocal since the 
vertices and propagators in Eqs. (4.12) and (4.22) are dependent not only on the relative 
three-momenta involved but also on their total three-momentum. They include the 
nonstatic (recoil) effects in all orders of the socalled ^ expansion . 



4-4 Other Clothed Eigenstates of Meson-Fermion System 

Up to now we have focused upon the clothed tiN and A^A^ states. Let us discuss other 
clothed states. 

If we start with the same "zeroth" approximation to the total Hamiltonian, our de- 
scription of clothed vriV and NN states will be very similar to that given for vrA^ and A^A^ 
states. Actually, it is the case where one has to deal with the charge-conjugated states. 
Here, we mean the nucleon-antinucleon conjugation. 

A different situation holds in the case of clothed fermion-antifermion 
and two-meson states. In fact, the operator K\ contains the interactions K^mr ^ 
NN). Therefore, superpositions of the vrvr states aJac^Q and the NN states bjdj^ taken 
separately cannot be Kza eigenvectors. So, one has to consider the eigenstates of a mixed 
kind, 

$ = y dkidk2(^mTiki;k2)aJ{ki)aJ{k2)^l 

+ J dp^dp2^NNiPi,ri;P2,r2)bJiPi,ri)dJ{P2,r2)n (4.25) 

Calculation of the scalar products of Kza^ = with {aJajQ \ and {bJdjQ \ leads 
to a set of coupled equations for the coefficients <l>,r7r and ^nn- Of course, one may obtain 
separate equations for each of them. Thereat, the eigenvalue equation for ^nn will involve 
some terms of the ^f'^-order. Obviously, to be consistent they should be disregarded within 
the ZA considered. 

In the analogous manner one can study the eigenvalue problem for clothed three- 
nucleon and more complicated states. However, handling with Kza, we enter into the 
3A^-problem only with the two-body interaction K{NN — > NN). It would be interesting 
to take into account the three-body (six-operator) interactions (irreducible to two-body 



26 



ones) that are present in the total Hamiltonian K starting from the term [i?, [i?, [i?, V]]] ~ 
g'^ not exphcitly written in Eq. (2.23). 

5. Possible Modifications of the Clothing Approach 

5.1 Heitler's Unitary Transformation 

The "clothing" in Sect. 2 has been realized in the framework of the Schrodinger picture, 
the "bare" and "clothed" operators being the Schrodinger ones. Heitler in his book ^ 
discussed the corresponding UT's in the framework of the interaction picture. To establish 
the relation with Heitler's approach we shall derive Heitler's equation (which determines 
his UT) starting from our equation (2.19) of Sect. 2: 

W{a)[HF{a) + Hi{a)]W\a) = K{a) = Ko{a) + Kj{a) (5.1) 

Remind that the free part KQ{a) of the total Hamiltonian K{a) (expressed in terms of 
the clothed operators a) is equal to Hpia), Hp^oi) being given by Eq. (2.6) in which bare 
operators a, h and d are replaced by the clothed ones Oc, &c and dc- Now, let us write Eq. 
(5.1) in terms of interaction picture operators defined, e.g., as 

ap(t) = e'^«(")*ape-'^«(")* (5.2) 

For this purpose multiply both parts of Eq. (5.1) by exp(zii'ot) from the left and by 
exp(— ifCot) from the right and use the notation 

A{t) = e'^'''A{a)e~'^°' = A(e*^o*ae-^-^«*) = A{a{t)) 

Then we obtain, 

W{t)[HF{t) + Hi{t)]W^{t) = Ko{t) + Kj{t) (5.3) 

Multiply both parts of this equation by W''{t) from the left and use W^W = 1 and 
Ko{t) = Hpif) (see above). We get 

[Hpit), W\t)] = ~Hj{t)W\t) + W\t)Ki{t) (5.4) 

The l.h.s. of this equation is equal to —i-^W'^{t) (see, e.g., Eq. (11.52) in [0). So, we 
find, 

Hi{t)W\t) - ^^W\t) = W\t)Kj{t) (5.5) 

This equation coincides with Heitler's equation written in |^ (Ch. 4, between Eqs. (15.6) 
and (15. 7i)). The relation of Heitler's notations and ours is 

S = W\ K = Ki, H = Hi 
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So, Heitler's basic equation is equivalent to our Eq. (2.19) or (5.1). But his goal differs 
from the goal of "clothing". He requires that Ki must not contain interaction terms 
which give rise to virtual processes. By definition, the latter can proceed in spite of the 
inequality of energies of the initial and final states (here "energy" means an eigenvalue 
of the free part of the total Hamiltonian) . Our bad terms also lead to virtual processes 
(e.g., VL —>■ 71 NN or N Ntt) but there are many other virtual processes not generated 
by bad terms (e.g., tttt NN or vrA^ —>■ ttttN at low initial energies). 

Of course, Heitler's requirement can be imposed also in the Schroedinger picture and 



this has been done by Sato et al. MB, |32 



5.2 Heitler-Sato Approach Versus the Clothing One 

We shall show here that under a condition all bad terms produce virtual processes. 
As there are many virtual processes which are not induced by bad terms, one may state 
that the Heitler-Sato condition is stronger than the bad terms elimination requirement. 

Let us remind the exact definition of bad terms. They are either two-, three- 
,. . . operator terms which contain only creation operators (and do not contain destruction 
operators) or three-, four-,. . . operators containing only one destruction operator. We call 
bad also the terms which are Hermitian conjugated to the above-mentioned. Using the 
notion of the class defined in Sect. 2 one may define the bad terms as terms of the class 
[n, 0] and [n, 1], n > 2 and their H.c. Note that two-operator terms of the kind a^{k)a{k) 
are not attached to the bad ones. 

Bad interaction terms are responsible for the processes 

no particles ^ 2,3,... particles 
one — particle 2, 3, . . . particles 

This property also may be considered as the bad terms definition. The processes (5.6) 
will be called bad below. 

Let us prove the Statement: All bad processes are virtual ones under a condition 
on the masses of interacting particles. Indeed, energy is evidently not conserved in the 
bad process "vacuum —>■ several particles": the initial energy is zero while the final state 
energy cannot be less than the sum of the final masses. Further, consider the bad process 
ttj ^ ai + 02 + . . ., in which the particle aj with the mass mj converts into particles with 
masses mi,m2, .... In the particle aj rest frame the initial energy is nij. The energy is 
not conserved trivially if sum of the masses Sjmj of the final particles exceeds nij. The 
energy is conserved if mj > T^imi and final particles possess nonzero momenta. The set 
of inequalities mj < T,imi, Wj is the very Statement condition on the particle masses. 
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In the Yukawa model the bad processes — A^vr, NNN, N Ntttc, etc. are 

certainly virtual, the process vr — NN being virtual if /i < 2m. Under this condition all 
the bad processes generated with the Yukawa model are virtual. 

Let us consider the simplest virtual processes of the Yukawa model which are of the 
order g^, namely Q —>■ NNtt, N —>■ Ntc, tt — > NN (if fi < 2m). They coincide with the 
(simplest) bad processes of the order g^. The corresponding three-operator interaction 
terms are of the order g^ and can be removed by the unitary operator W = exp R where 
i? is a three-operator expression of the order g^ (see Sect. 2). 

6 The Okubo Blockdiagonalization Method 

Here, following we regard the UT H — > Hu = U^HU that makes the Hamiltonian 
H (generally speaking, an Hermitian operator) block diagonal (cf. our brief discussion in 
Introduction). After the transformation the primary H eigenvalue problem is reduced to 
the diagonalization of separate neardiagonal blocks of Hu. 

It is well known that the transformation H Hjj can be interpreted either as the 
connection between the two matrices of one and the same operator H with respect to the 
different orthonormal bases, or as the relationship between the two operators H and Hjj 
unitarily connected via the operator U. We shall start from the first point of view. 

6.1 UT as Change of Basis 

Let Hnin = {n'\H\n) {\/n',n) be the H matrix with respect to a complete set of or- 
thonormal vectors \n) and H,y/i, = (z/'|if|z/) (Vi/', z/) represents H in another orthonormal 
basis |z/). The indices n{i>) can take on discrete or/and continuous values. One has (see, 
e.g., |§, Ch.I) 

if^v = (i^'lHlu) = Sn'Sn {iy'\n'){n'\H\n){n\iy) (6.1) 

where S„ denotes a sum or/and an integral over n. 

The r.h.s. of Eq. (6.1) can be written as the matrix product 

S„.S„ {W).'n'Hn'nUnu = {WHU),>, (6.2) 

where we have introduced the notation Unu = Then 

(f/t),,„. = u:,,, = {nyy = [vY) 

The transition matrix U is unitary in the following sense: 

(U'^U)^'^ = S„ {u'\n){n\u) = 5^>y (6.3a) 
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and 

{UW)n'n = {n\u){u\n) = {n'\n) = 6n'n (6.36) 

Recall that all these basis vectors are orthonormal. 

Whenever U is given one can express |z/) in terms of |n), 

|z/) =SJn)f/„, (6.4) 

One should point out that the spectrum of indices n enumerating the vectors \n) 
needs not to be identical to the spectrum of indices u enumerating |z/). For example, let 
us consider the Hamiltonian Hos of one-dimensional quantum oscillator, whose matrix 
in the coordinate representation can be made diagonal by using the Hermite functions 
M^) (z/ = 0,1,2,...) : 

oo oo 

{iy'\Hos\i^) = / dx' / dx^p*,{x') {x'\Hos\x)^pi,{x) ^ 6^1^ (6.5) 



(cf. Eqs. (6.1)-(6.2)). In this case the transition matrix (xlu) = ipv{x) is not "square" 
because its columns(rows) are enumerated by the integer(continuous) numbers. Such a 
matrix can be called rectangular. Note that in a finite-dimensional space the transition 
matrix always is square. This exemplifies that one cannot, in general, diagonalize H by 
means of a square matrix. 

In this connection, one should note that the diagonalization scheme considered by 
Okubo and exposed below is not aimed at a perfect solution of if-eigenvalue problem. 
Rather, it sets a more humble task, viz., to find (approximately) some of its eigenvalues 
and eigenvectors. Okubo suggested a realization of the scheme via a square matrix U . 
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6.2 Blockdiagonalization in Matrix Form 



In accordance with [0 we require that the H matrix in the new basis 

H,.,= {U^ HU),,, (6.6) 
should have the blockdiagonal form, 

U^HU = K=\ " I (6.7) 




with the two blocks Kn and K22 the meaning of which will be clarified a little later. This 
requirement must determine the matrix U that in its turn enables one to construct the 
new basis vectors |i/) (see Eq. (6.4)). 

Now, following Okubo, we shall confine ourselves to finding such a square matrix 
U that \v) = U\n). In other words , it is assumed that one can set an one-to-one 
correspondence between the indices n and u. 

First of all, let us turn to the definition of Kn. It is a block of the matrix K with 
elements that are enumerated by indices i>i (or ni) which belong to a subset of all u (or 
n) values. The K22 elements are enumerated by the remaining u values denoted through 
h'2- Let us give examples of the ui choice. 

One may take one value of z/ as z/i, viz., the index of the vacuum state. Hence, Ku 
has one element. If we are able to find U which leads to Eq. (6.7) then we can construct a 
normalized H eigenstate, namely the physical vacuum |0) = S„ |n) f/„o- After this step we 
consider the block 7^22 as a starting matrix for the subsequent blockdiagonalization, viz., 
to introduce a new set of z/^ (e.g., let u'^ be indices of the one-particle Hp eigenstates), to 
find a new U'. At the next stage one can enumerate elements of the recurrent block 11 by 
indices of two-particle states, e.g., the states "two nucleons, no mesons" (this is Okubo's 
example, see Sect. 2 in |^). One more choice of Ku will be discussed in Subsect. 6.5. 

The option of z/i or ni allows one to divide H into the four blocks: Hu with elements 
{n[\H\ni), H22 (elements {n2\H\n2)), H12 (elements {ni\H\n2) ), and H21 = Hi2^ . The 
matrix U can be represented analogously. Keeping this in mind the l.h.s. of Eq. (6.7) 
can be rewritten as the product of matrices composed of the blocks described above (see. 



e.g., Ch. 0.7 in and Ch. 1.6 in [P 



Ul2^ U22^ ) \ ri2\ -H22 / V ^21 U22 

One can consider Eq. (6.7' ) as the equation for U. Okubo suggested to seek its 
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solution in the class of unitary matrices of the kind (see Note underneath) 

\U2i U22 ) \A{h + A^A)~-2 {h + AA^)~-2 j 
where /i(/2) is the unit matrix for the subset of indices ^1(712), and A is an rectangular 
matrix with elements A^jni of the U21 kind, which should be determined. Under this 
convention the unit matrix for the full set of indices has the block structure 

1=1 ' 1 , (6.9) 



/. 



2 



and we have the properties 

All = l2A = A , (6.10a) 
[A^A,Ii]=0 (6.106) 
The matrices A'^A (of the Un kind) and AA^ (of the U22 kind) are square, hermitian 



and positively definite, so that the square roots y h + A^A and \/l2 + AA^ can be defined 



(see, e.g., Ch. 5.8 in |^). It is readily verified that the r.h.s. of (6.8) meets left unitarity 
U"^ U = 1. Other unitarity condition UW = 1 can be proved by means of equality 
Af{A^A) = f{AA^)A that is valid if f{x) is a polynomial or a series of x. 

The requirement for blockdiagonalization, K21 = (or K12 = 0), gives the equation 

H21 + H22A - AHu - AH12A = (6.11) 

for A determination. This equation is equivalent to the condition (13) from [|^. 

Eq. (6.11) is nonlinear and it can be solved exactly only in a few simple cases (see 
Sect. 6 in and Appendix C ). For realistic field models one has to develop a perturbative 
method in order to find A (see Subsect. 6.5). 

Note . Generally speaking, the left unitarity for a matrix U with the nonsingular block 
Uii enables to express blocks Uij{i,j = 1,2) through the matrix A such that U21 = AUu 
and blocks = Si and 5*22 = S2 of an arbitrary matrix S, 

Si 

^2 

with Si^Si = Ii and S'2^S'2 = I2 (cf. Eqs. (10) in 0). The corresponding matrix Uokubo 
(i.e., its representation after Okubo) can be written as the product US = Uokubo where 
U is given by (6.8). However, since the UT via 5* conserves the blockdiagonal structure 
of K it is sufficient to consider only the form (6.8). Of course, Okubo's representation is 
not most general form for unitary matrices. The simplest exception is the 2x2 matrix 

1 

1 



S 
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6.3 Blockdiagonalization with Projection Operators 



Here the UT H' = U"^ H U is interpreted as a relation between the different operators 
H and Hf^, where ?7 is a unitary operator. By definition, the operator is a hnear one- 
to-one mapping of a Hilbert space Ti onto itself, i.e., an isomorphism of Ti. In particular, 
it transforms the orthonormal basis vectors \n) ^Ti into the vectors \n) = U\n) of other 
orthonormal basis G Ti. Therefore, we can write 

{n'\H'\n) = {n'\WHU\n) = {n'\H\n) , (6.12) 

Eq. (6.12) sets up close links of H' = U'^HU with the UT considered in Subsects. 6.1 and 
6.2. 

One can deal with H' = U^HU in the same manner as in Subsect. 6.2 introducine; 
matrices of H, U, H'^ with respect to a basis \n) and representing them in block form, eq. 
(6.7 '). But one can realize the Hamiltonian blockdiagonalization without reference to any 
basis on Ti. We shall introduce (cf. |0) the projection operators rji and rj2 = 1 — rji onto 
a subspace Hi C H and its complement H2 such that H = 7-^i07-^2H- Mathematically 
strict definitions of projection operators and their properties in a Hilbert space can be 
found in (Ch.2). Of great importance for us is the property 

ViVk = VkVi = ^ikVi (6-13) 

(i,k = 1,2) 

Then, for any operator O in 7i one can write 

= {vi+V2)0{7]i + r]2)= E O*'' (6-14) 

where 

0'^ = riiOrij (6.15) 

with 

rj, O"' = 6,kO'\ O"^ VI = 0''6,i (6.16) 

In terms of such decompositions the product of the two operators A = J2ij and 
B = -B*"' can be written as 

AB = Y,A'^B^'' , (6.17) 

ijk 



-"^^Unlike the UT K — W HW^ in Subsect. 2.3 that represents the same operator H. 

^^Of course, the particular option of Tii depends on the nature of the problem (see Appendix C and 

Subsect. 6.5). 
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so that 

{ABY'' = Y,A"^B^'' (6.18) 
j 

By using the rule we find 

H' = WHU = {U^)"''W^U^^ (6.19a) 

imjk 

or 

H' =J2 {W^'YH'^W' , (6.196) 

imjk 

Remind that the projectors rji are hermitian , i.e., rji^ = rji. 

One can set an one-to-one correspondence between the operators if*-' and the matrix 
blocks Hij introduced in Subsect. 6.2. Actually, each operator H^^ = 1,2) acts in 
the full Hilbert space Ti. but its matrix with respect to the basis {\n)} has merely one 
nonzero block Hij with the three remaining blocks being zero. Analogous relation takes 
place between and Uij. 

Further, the operators in Eqs. (6.19) can be expressed through an arbitrary 
operator A of the kind 21, i.e., A = f]2Ar]i. For instance, one can write (cf. Eq. (6.8)) 

[/2i = ^f/ii = ^(1 + AU)-^7]i = A{1 + AU)-^ , (6.20) 

Okubo's requirement for decoupling the two subspaces Hi and 0.2 means that it must 
be H'^^ = or according to Eq. (6.19a) 

^(f/t)2i^iif/ii = (6.21) 

It leads to the nonlinear equation for A 

r]2{H + [H, A] - AHA}r]i = (6.22) 

As a matter of fact, it is the same equation as Eq. (6.11), viz., one may consider the 
latter as the record of Eq. (6.22) for the only nonzero block A of the operator A. 

The solution of Eq. (6.22) , if it exists (cf. the discussion in ), yields the Hermitian 
operator 

H'^^ = r]i{l+A^A)-^l + A^)H{l + A)il+A^A)---r]i , (6.23) 

The resultant operator H' given by Eqs. (6.19) (see also Eq. (6.23)) has a block diag- 
onal structure that simplifies determination of its eigenvalues and eigenvectors. However, 
H' H and therefore H' eigenvectors are not H eigenvectors. Recall that our UT is aimed 
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at solving the H eigenvector problem (at least, its partial solution). But if we have suc- 
ceeded in finding H'^^ eigenvectors one can easily get the corresponding H eigenvectors. 
Actually, if we have H'^^^>i = E^^i G Tii), then is also H' eigenvector, 

H'^^i = (i/'" + = (6.24) 

It follows from U'^ HU^i = E^^i that f/\l/i is H eigenvector, 

HU'^i = EU^i (6.25) 

The obtained relation between H' and H eigenvectors can be represented in another 
form. Let us expand the H' eigenvector in the basis vectors |?7,), 

l^l) = Sn\n)Cn 

The corresponding H eigenvector f/^i has the same expansion coefficients c„ with respect 
to the other basis \n) = U\n), 

f/^l = SnU\n)cn = Sn\n)cn (6.26) 

Exact solutions of the decoupling equation (6.22) can be derived for a few simple cases 
(see, e.g.. Sect. 6 in 0|). One of them is the so-called scalar field model (see Ch. 12 in 
||15|| ). The respective solution is very instructive having many attributes of more realistic 
field models. It is given in Appendix C. 



6.4 New Creation-Destruction Operators within the UT Method. 
Comparison with the Clothing Procedure 

As before (see Sect. 2) we consider that the original field Hamiltonian if is a function 
(functional) of the bare creation-destruction operators. We denote their set by the same 
symbol a, ap being one of them. Moreover, one may assume that the unitary operator U 
also is a function of a. Then, H' = U'^HU may be written as 

H'{a) = W{a)H{a)U{a) (6.27) 

In order to compare the Okubo approach and the clothing procedure developed in 
Sect. 2 let us introduce the set a of new creation-destruction operators (with the denota- 
tion ap for one of them) defined as 

ttp = V{a)apV^{a) \/p , (6.28) 
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or 

ap = V^{a)dpV{a) , (6.29) 

where V{a) — V{d) 7^ 1 is an arbitrary unitary operator (not necessarily coincident with 
U{a)). 

Applying the UT V to the both sides of Eq. (6.27), we have 

V{a) H'{a) V\a) = H'{a) = U\a) H{a) U{a) (6.30) 

Further, by means of (6.29) the total Hamiltonian can be expressed (cf. the derivation 
of Eq.(2.17) ) in terms of the new operators a, 

H{a) = H{V\a) d V{d)) = V\d) H{d) V{d) (6.31) 

One should stress that the operators H'{d) and H{d) are different from H'{a) and H{a), 
respectively, \IV ^1. However, \IV — U, than H'{d) turns out to be equal to the starting 
Hamiltonian H{a). Actually, in this case the r.h.s of Eq. (6.31) gets to be equal to H'{d) 
in accordance with Eq. (6.30), so that 

H{a) = H'{d) = U\d)H{d)U{d) (6.31') 

The option V = U = gives rise to the clothed operators a = a = W^aW (see 
Subsect. 2.3), and then H\a) = K{a) = WH{a)Wl 

Thus, our consideration shows how the UT H ^ H' = U^HU can be reduced to a 
transformation of clothing type. Of course, there are distinctions between the Okubo 
approach and the procedure shown in Sect. 2. They are due to those purposes which are 
inherent to each of them. We shall return to this point in Subsect. 6.6. However, let 
us note here that although the operator Uokubo{ci) determined by solving the decoupling 
equation (6.22) and the operator U ciothing{o) — W^a) arc different functions of a one can- 
not d priori exclude a resemblance or even the perfect coincidence of some approximations 
to them. 

Note one more definition of d by the relations 

{n'\dp\n) = {n'\ap\n) \/n',n,p , (6.32) 

where |n) = U \n) (cf. the beginning of Subsect. 6.3). In particular, Eq.(6.32) means that 
the matrix elements of the new meson destruction operator d{k) with respect to the new 

vectors |n) are equal to the corresponding matrix elements of the bare meson destruction 
operator a{k) with respect to the old vectors \n). It is clear that this definition with 
U = U {a) is equivalent to (6.28). 
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6.5 Perturbative Construction of Okubo's Unitary Transformation. 
Elimination of Mesonic Degrees of Freedom 



Here we present the approximate solution of Okubo's equation and determination of 
the block H'^^ that have been given in [^. Just as in that paper let us consider the 
system of fermions (nucleons) and mesons (pions) with Yukawa-type interaction linear in 
the meson field (see Sect. 2). In this model the basis states \n) are enumerated by the 
two indices \n) = \mf), viz., m{f) enumerates the meson (fermion) states, \Omf) being 
the state without mesons (the NN contents may be arbitrary). In accordance with the 
Okubo idea we decompose the full space H of meson-nucleon states into two subspaces 
(sectors), namely, the fermion (nucleon) sector Tio which is composed of no-meson states, 
being spanned onto the subset | 0^/), and its orthogonal complement Ticomp that consists 
of the states with nonzero meson number. The projection operator rji into T^o ^^"^ t)e 
constructed as 

r/i = S/ |0„/)(0„/| . (6.33) 

By definition, 

«(fc)|0^/) = V/,fc, (6.34) 

and therefore 

a(fc)r/i = T]ia\k) = (6.35) 

The operator of interest H'^^ has the following matrix structure (cf. the note after 
Eqs. (6.19)): 

^ ( : ) . (e.3e) 

where the block H^-^^ consists of the elements {Omf \ H' \ Omf), V/', /. At the same time 
the matrix 

H"' = I ° ° I , (6.37) 

contains the elements {m'f \ H' \ mf), Vm' ^ Om ^ rn , V/', /. For brevity , henceforth 
the index 0^ of the no-meson state will be replaced by 0. 

In order to simplify the subsequent equations let us separate the matrix of any operator 
O that acts in 7i = Hq © 'Hcomp into the subblocks [Om/m] such that 

[Om'm]f,f = {m'f \0\mf) (6.38) 
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In pa the letter P has been used instead of t]i. 
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So, the elements of these subblocks are marked by the fermionic indices /. The subblocks 
can be called fermionic. The block Hu coincides with [Hqo] , etc. The rectangular matrix 
A, which determines U in Eq. (6.8), consists of the subblocks [Am/o] with m' 7^ 0, i.e., 
the submatrices Aj^^ with m' = k for one-meson states marked by the momentum k , 
the submatrices [^fc^j^^ o] "^i^^ — (^1^2) for two-meson states , and so on. 

Now, the key equation H' = U'^HU can be rewritten as the set of equations for the 

subblocks [f/m'm] 5 

[H'^'j} = Srn'Sm[U\irn'][Hjn'rn\\Um.fM]-, V/i', /i , (6.39) 

where runs the same values as m does. 

In their turn the subblocks in question may be considered as the matrices of operators 
which act onto fermionic degrees of freedom. Such operators have been represented in 



25| as functions of the bare creation-destruction operators but fermion ones only. Let us 



denote the set of them by the symbol /. Then, one can write, e.g., 

[U„,'mCf)]f'f = {m'f\U\mf) (6.40) 
Similarly we introduce the operators H'^i^{f) {Hm'mif)) which correspond to 

Under such convention the relation H' = U'^HU can be reduced to a hybrid form 
being expressed in terms of the functions of / and the matrix elements with respect to 
the meson states, 

= S^^S^U\,^\f)H^,M)Un.,{f) (6.41) 

Such representation of the operators turns out to be convenient when finding approx- 
imate solutions of Eq. (6.22) which can be rewritten as 

r]2{Hi + [Hp, A] + [Hi, A] - AHiA}r]^ = (6.42) 

since rj2HFi]i = 0. 

We consider its approximate solution assuming 

A = + A^'\A^'^ ~ (7°, A^'^ ~ (6.43) 

The contribution A^^^ is absent since with interaction switched off (g 0) the Hamilto- 
nian H = Hp is already blockdiagonal. In this connection, remind that Hj = V + Mren 
where the mass counterterms Mren = Mferm + Mmes are determined by Eq. (2.7). 

Such A = A^^'^ turns into zero only those terms of Eq. (6.42) which are of the ^f^-order, 

i.e., 

r]2{[Hp,A] + V)r], = , (6.44) 
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One can consider Eq. (6.44) as a relaxed form of Okubo's constraint H'^^ = (or Eq. 
(6.42) ), which may be imposed instead of the latter. 

Using the commutativity of Hp and rji let us rewrite Eq. (6.44) as 



(6.45) 



This equation is of the same type as Eq. (2.21) for the operator R, viz., [Hp, R] —V — 0, 
whose solution is given in App.A. Therefore, taking into account formula (A. 4) we find. 



A^ilim^j V^^{t) exp {-€t)dt 



(6.46) 



For interaction V linear in the meson field we have Ajno{f) = if the index m corre- 
sponds to two- , three- , . . . meson states. Let us consider the nonzero subblock Aj^^ , 
i.e., the subblock with elements {kf \ A \ Of). With the help of the representation (A. 9) 
for the interaction V and relation 



{a^{k)n I a^{k',t) \ Q) = 5{k - k') expicoj^t 



we obtain 



with 



^koif)-\^xI^kom) 





exp 



(o;^ + te)t dt , 



^fco(/W) = FKt)V-f'F{t) , 



(6.47) 



(6.48) 



where f{t) — exp iHFfgrmif)^ / exp —iHpferm{f)t is the subset of the fermion opera- 
tors in the interaction picture. We imply the division Hp — Hprnesi'Th) + Hpfermif) into 
the mesonic and nucleonic parts (see Eq.(2.8)), m being the subset of the meson operators 
a{k) and a^{k) Vfe. For notations F and see App. A. 

At the same time the (fcO) subblock of the operator R — TZ — T?) determined by Eqs. 
(A.4) and (A.9) is 



dt 



(6.49) 



Comparing Eqs. (6.47) and (6.50), we get the solution for ^4^^, 



(6.50) 



in a compact form. 

Here, we do not intend to write down explicit expressions for these subblocks in terms 
of the fermion creation(destruction) operators. The respective result for A^^^[f) with 
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some extension to other Yukawa-type meson-nucleon couplings can be found in (see 
Eqs.(17)-(18) therein). 

After all, in order to derive an approximation to the operator H'^^, let us employ the 
expansion of (1 + A^A)~^ in the powers of A'^A ~ (?^, 

^ 1- i^U (6.51) 

Then, we obtain from Eq. (6.23) the expression 

H"^ = r]i{H + A^H + HA+ ^A^[H,A\- ^[H,A^A}r]i + 0{g^), (6.52) 

from which and condition (6.22) for A it follows that 

H'^' = Hpf,Uf)Vi + mIIM)Vi + lvi{A^H + HA}vi + 0{g') , (6.53) 
or neglecting the terms of order g^, 

H"' = Hpf,rmif)Vi + MflXf)Vi + \vi{A^'^^V + V^(i)}r7i (6.54) 
(cf. Eq.(15) of |]25|). When deriving Eq. (6.53) we have used the equation 

r^,Hr^, = Hpf,rmif)m + MgL(/)r/i + Oig') , 
whereas the change from (6.53) to (6.54) is based upon the properties 

r]2Hr]i = ri2M'^l^{m)rii + 772^^ + 0{g^) 
r]iHr}2 = riiM'^l^{rh)ri2 + r/iKr^a + 0{g^) 

and the fact that 

ilrM^2s{^)Avi = 0{g') , 
r^^VAiii = r^iVA^^\^ + 0{g^) , 

for V linear in the meson field. 

The corresponding subblock can be expressed through the subblocks Ap,^ and 

H[, = HpfeUf) + Mfl^if) + \j{A^ okUWkoU) + H.c.}dk = 

HpfermCf) + Mfl^) + \j {F^ F ■ F^'^ F + H.c.}dk (6.55) 
where we have employed Eq. (6.50) and Eq. (6.48) at t = and their H.c. 
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6.6 Clothing Procedure vs Okubo Approach 



Let us discuss in detail common and distinctive features of the two kinds of UT. Partly, 
we have concerned this subject in Subsect. 6.4 ( see also App. C ). 

The Okubo and clothing approaches differ in their goals. The UT by Okubo is aimed 
at to nullify the nondiagonal blocks ifgi and of the transformed Hamiltonian. The 
clothing and akin approaches (see Sect. 5) require that some undesirable operator terms 
("bad" or "virtual") must be absent in the transformed Hamiltonian. 

In the context it is worthwile to mention Nishijima's modification of the Okubo 
idea. Instead of Eq. (6.8) he used the representation U = . . . exp {182) exp (iSi), similar to 
that employed within the clothing approach. But hermitian operators Sn {n = 1, 2, 3, . . .) 
are to implement Okubo's requirement, viz., the operator Sn must remove K21 matrix 
elements of the ^f^-order (instead of removing some bad terms ~ g""). 

No wonder different approaches give, in general, different resulting transformed Hamil- 
tonians. Nevertheless, we want to show that the both approaches can give some coincident 
results (even for the realistic field model), being realized approximately (in a perturbative 
way). In particular, effective quasipotentials of the g^ order for nucleon-nucleon interac- 
tion turn out to be the same. 

Before comparison of the resulting Hamiltonians we must remind that Okubo's H' 
(see Subsect 6.5) does not coincide with the starting Hamiltonian H whereas the operator 
K of the clothing approach does. H' can be considered as a function of bare operators 
a while K depends on clothed ones a. However, we have shown in Subsect. 6.4 that H 
coincides with the operator H'{a) which depends upon new destruction-creation operators 
a = UaW in the same manner as H' depends on a. So, H'{d) and K represent the same 
operator and can be compared to each other. Let us recall, however, that a may not 
coincide with clothed operators a, and H' and K may not be the same function of their 
arguments. In what follows we imply that H' means H'{a) and is a function of 
fermion destruction-creation operators / = UfUK 

Now, we can prove the statement: given by Eq. (6.55) coincides with the fermionic 
part of Kp + M^ll + V] (see Eq. (2.23)), i.e., with 

Kferm = Kpferm + M)J^^ + (fermiouic part of -[i?, V]) (6.56) 

More exactly, H'^^ and Kferm are the same function of their arguments (fermionic 
operators / and fc, respectively). Indeed, the last term in Eq. (6.56) is that part of 
V] ~ g^ which depends only on fermionic operators. We calculate V] in App. A. 
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Its fermionic part (denoted below as FP) is contained in |[7^, V^] + H.c, for [7^, V^] see 
the last term in the r.h.s. of Eq. (A. 14). We have 

FP=^J [F^R^F ■ F^V-^F + H.c.]dk (6.57) 

Comparing this expression with the r.h.s. of Eq. (6.55), we arrive to the above statement. 

The fermion mass counterterm Mjg|.^ in Eqs. (6.55) or (6.56) must cancel two-operator 
fermionic contributions which arise after normal ordering of FP (see Subsect. 2.5). One 
should note that there are no other such terms in y]. So, we may omit Mj^^.^ from 
(6.55) or (6.56) provided FP is replaced by the normally ordered counterpart : FP : of 
(6.57). 

Let us mention that similar evaluation of M^l^ (cf. Subsect. (2.5)) would take addi- 
tional efforts within Okubo's approach (see Note at the end of this subsection). 

Of course, we may explicitly express FP given by Eq. (6.57) through b, d, and d^ 
opening the abbreviations accepted for F'^V'^F and F'^R^R in App. A (we omit here the 
symbol tilde or the subscript c when handling with the fermion operators). Such repre- 
sentation shows that : FP : contains "bad" terms of the kind fe^rf^fe^c?^, b'^d'^b'^b and b^d'^d'^d. 
As stressed in Subsect. 2.4., they should be eliminated from K via the clothing transfor- 
mation W4. Such terms are also unpleasant within Okubo's approach, viz., they prevent 
the no-fermion state Qf and onc-fcrmion ones b^flf and d^flf to be H'n eigenvectors along 
with two- fermion states of the kind (4.14). 

In other words, the problem of finding H[i eigenvectors is not essentially easier than 
that for the starting Hamiltonian. In the spirit of Okubo's approach the elimination of 
"bad" terms can be implemented by performing UT of H[i such that the transformed 
subhamiltonian (i^d)' = U'^H[^U' would not contain matrix elements corresponding to 
the processes Jl/ — > two pairs, N — > A'"-|-pair, etc.. The relevant projector ri[ may project 
on states without pairs. 

As argued in Sect. 4, these additional transformations W4 or U' do not alter the 
remaining "good" four-fermionic pieces of : FP : which are of the g"^ order (see, e.g., Eq. 
(4.18)). 

So we obtain from H[i or Kferm the same expression Hpferm + G : FP :, where 
G : FP : is a "good" part of : FP :. Using the notations of Sect. 4 (see Eqs. (4.1)-(4.3)), 
it may be written as (Kza) ferm or 

K2ferm + K{NN N N) + K{NN NN) + K{NN NN) (6.58) 

The no-fermion, onc-fermion and two-fcrmion deuteron-likc eigenvectors of this operator 
give some approximations to the corresponding H eigenvectors. 
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Note . The zeroth approximation operator Kza considered in Sect. 4 enables us to 
find also approximations to H eigenvectors which describe one-meson states and meson- 
nucleon bound and scattering states. In the Okubo-Korchin-Shebeko approach exposed 
in Subsect. 6.5 one should use for this purpose the submatrix H22- However, the latter 
embodies, in general, undesirable matrix elements which correspond to the processes one- 
meson other states (e.g., two meson + pair), etc. In order to simplify the finding of 
H eigenvectors under discussion one needs additional Okubo's transformation of the kind 
{H'22)' = WH22U such that {H'22)' would have no undesirable matrix elements. 



7 Conclusions 



A considerable part of this work is devoted to development of the UT method using 
the so-called clothing procedure in RQFT. This procedure has two aspects. On the one 
hand, we express with its aid a total Hamiltonian H for interacting fields in terms of the 
new operators which correspond to the creation (destruction) of clothed particles. The 
latter possess, by definition, the properties of observed (physical) particles and are as 
a matter of fact " quasiparticles" within our approach, if one draws a parallel with the 
method of canonical transformations in quantum theory (see, e.g., §§52, 84 ). On the 
other hand, such representation for H, being implemented partially or perfectly, enables 
us to formulate an approach to solution of the H eigenvalue problem. 

We have discussed in detail various UT of H in Sect. 6. In this connection, we 
distinguish two kinds of UT's destined for approximate determination of H eigenstates. 
Both may be deduced from the relation between the Hamiltonian H matrices with respect 
to two different sets of basis vectors . 

The first kind may be written as H' = HU, where H is the input Hamiltonian 
which is subject to UT, the transformed operator H' being not equal to the input one. 
However, H eigenstates can be obtained from those of H' using U (see Subsect. 6.3). 

The second kind is determined as H = u'^ HU {U = w'^). Now, the transformed 
Hamiltonian is the input one, but differently represented, while the operator H ^ H is 
related to if in a simple way. These kinds of UT's are used in literature but authors 
sometimes overlook that H' eigenstates do not coincide with those of H. 

Applying either of such UT's one may impose different constraints on the form of 
the transformed Hamiltonian that leads to the definite recipes for constructing the cor- 
responding unitary operators. We have shown the resemblances and distinctions of some 
known applications of the UT method: the clothing procedure, the approach by Heitler 
and Sato, and Okubo's blockdiagonalization (Sects. 5 and 6). One should note that in 
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all applications the unitary operators U is determined approximately (with the exception 
of simple models). . 

Our consideration of the problem of bound and scattering states differs from akin 
approaches (cf., 0, ^5[) with the following distinctive feature. Our interactions between 
clothed particles ( as an illustration, the vrA^ and A^A^ quasipotentials) are parts of a single 
operator Kza which can be regarded as a zeroth approximation to the total Hamiltonian. 
If the A^A^ system is considered then Kza generates the A^A^ Hamiltonian, whereas for 
the ttN system Kza gives rise to the vrN Hamiltonian. This makes clear the relation of 
the eigenstates of such Hamiltonian to the eigenstates of the original field Hamiltonian: 
the former are some approximations to the latter. 

One should point out that nonlocal properties of these quasipotentials have a relativis- 
tic origin. In this context, we would like to note a growing interest in similar relativistic 
effects in the modern theories of nuclear forces(see, e.g., |56|, |57| and refs. therein ). 
We show how the mass renormalization program is realized within the approach devel- 
oped here. The transformed Hamiltonian when expressed in terms of clothed operators 
turns out to be dependent on renormalized (physical) masses and not bare ones. Some 
tricks shown in App. A can be useful in future calculations of the radiative corrections 
(renormalizations ) to bare (trial) masses for field models with a cutoff in the momentum 
representation. 

We prove that clothing UT's remove the undesired ("bad") terms simultaneously from 
the Hamiltonian and the generators of Lorentz boosts. Our proof is valid for any RQFT 
model (see Sect. 3). 

Our three-dimensional formalism is covariant in that sense that we give definite pre- 
scriptions for the transformation properties of the clothed states with respect to the 
Lorentz boosts. 

With the help of a simple example, we demonstrate in App. B that the clothing 
transformation W of Sect. 2 may happen to be implemented not by an unitary operator 
in its usual mathematical sense. We argue that W ought to be understood as an element 
of some algebra (lacking any operator representation), being unitary in an algebraic sense. 

At last, we show in App. C original tricks and results concerning a nonperturbative 
solution of Okubo's decoupling equation. 

We believe that the concept of clothed particles and the approach exposed here can 
be applied to different areas of the nuclear physics: from the theory of nuclear structure 
to description of nuclear reactions including the processes with meson production. 
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Appendix A 



Three-operator clothing transformation for the Yukawa model. 
Four-operator interactions between clothed particles and their 

normal ordering. 

The defining equations for the Yukawa model are given in Sect. 2. We use throughout 
this article notations of [^] assuming, in particular, = •yolfilo (a* = 0,1,2,3), q = 



q^l^l, 75 = i7°7W^7^ = 75 and 7I = 1. 

A.l. Three-operator clothing transformation W = expR, see Sect. 2.4, can be 
found by solving the equation [R, Hp] + V = 0. The antihermitian operator R is assumed 
to be of the form R = TZ — TZ'^ with TZ given by Eq. (2.22). The commutator [R,Hp] 
can be directly evaluated by using the commutation relations (2.5). Then we obtain the 
equations for the coefficients Rij involved in TZ, see Eq. (2.22). Their solutions are 



Rii{p'r';pr) 



vtdp'r'-pr)l{Ep,-Ep 



R%{p'r';pr) - V,^^{p'r';pr)/iEj^, + Ep - uj^) , 
R^,{p'r';pr) - vi{p'r'-pr)/{-Ep, - Ep - ujf,) 
RUp'r';pr) = vM{p'r';pr)/{-Ep, + Ep - u:^) 

Here we have used the notations 

V^^{p'r']pr) V^{p'r']pr)^ ^ g m 
V^^{p'r'-pr) V^^{p'r'-pr)\~'WW'^j2^^E^ 

u{p',r')'j5u{p,r) u{p',r')'j5v{-p,r) 
_ v{—p', r')-fr,u{p, r) v{—p', r')-f5v{—p, r) _ 



(A.l) 



X 



6{p + k- p') 



V^{p'r';pr) . 



(A.2) 



Eqs. (A.l) have meaning if the denominators in their r.h.s. are not zero (note that 
p-\-k — p' = according to (A.2)). One can show that this is the case under the condition 



fi < 2m . 



(A.S) 



The physical sense of this condition is discussed in Subsect. 2.4. 
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Alternatively, the solution R of Eq. (2.21) can be represented as 

/•oo 

R=-i lim / V{t)e-''dt , (A.4) 

where V(t) = exp{iHpt)Vexp{—iHpt) is the interaction operator in the Dirac picture. 
Obviously, V{t) is determined by Eq. (2.16) where the operators ac{k), bc{p,r) and 
dc{p,r) are replaced by ac(fc)exp(-icj^t), bc{p,r)exp{-iEpt) and dcip,r)ex-p{-iEpt), respec- 
tively. Therefore, the evaluation of R is reduced to integrals of the kind 

/ e*("+*^)* dt = ^— , (A.5) 

JO X + ie 

where x is any of the denominators in (A.l) and e > 0. The limit e — > 0+ in Eq. (A.4) 
exists, and it is finite if the inequality (A. 3) takes place, i.e., if x 7^ 0. This evaluation of 
R shows readily that the solution given by Eq. (A.4) coincides with that determined by 
Eqs. (A.l). Also, one can directly verify that (A.4) meets Eq. (2.21). In fact, we have 

roo 

[R, Hp] = -i ^hm e""* [V{t) , Hp] dt 

= lim / e""*— ^(^) dt = -V . (A.6) 

The last equality follows (under the condition /i < 2m) if one calculates ^V{t) and then 
integrates with the help of (A.5). We use the form (A.4) in Sect. 3 when constructing 
the Lorentz boosts in terms of the clothed operators. 

A. 2. Evaluation of [i?, V] is a tedious exercise that can be simplified with the aid 
of a more compact notation. Indeed, expressions (2.16) and (2.22) for V and R have the 
identical operator structure, viz., 



blip', r')xfi(pV'; pr)6c(p, r) + b^p', r')X^2{p'r'- pr)ai'-^^- [-p, r) 

+ d,{-p',r')X^^{pV;pr)b,ip,r) + d,{-p',r')X^2{p'r'-pr)dl{-p,r) (A.7) 
with definite c-number coefficients xj^ = 1,2). 



Now, let us rewrite (A.7) matrix product 

" Xf^ (p V; pr) X^2 {p'r'] pr) 
. X^i {P'r'] pr) X^2 iP'r'; pr) 



bciP,r) 
dl{-p,r 



{bl{p',r'), d,{-p',r')) 

= F^{p',r')X^{pV;pr)F{p,r) , {Al 



According to the prescription given by Eq.(A.4) the positive parameter e should be put equal to zero 
at the end of the calculations. 
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where along with the 2x2 matrix composed of the coefficients X^^ (cf. Eq. (A. 2)) we 
have introduced the operator column F and row 

nP.r) = , F\p',r') = iblip',r'), 4(-p',r')) ■ 

\dli-p, r)J 

The subsequent operations become even more concise after adopting the convention 

j dp' j dp Y.F\P' ^^')^^^P'^''^P^)P^P^^) = F^X^F ■ (A.8') 

r'r 

Under these notations, Eqs. (2.16) and (2.22) look as 

Jdk fW'F [a^{k) + al{-k)] 

^ Jdk F^V^Fa^{k) + H.c. = V + , (A.9) 

R = Jdk F^R^Fa^{k) - H.c. = TZ-'R) , 
where H.c. means taking the Hermitian conjugate of the first terms. Note that (y'^)t = 

A. 3. After these preliminaries, we have 

[R, v]^[n- n\ V] = [7^, v] + [v\ 7^t] = [7^, v] + h.c. (a.io) 

Keeping this in mind, it is sufficient to evaluate 

[7^, = [7^, V] + [7^, vt] . (A.ii) 

By using the definitions from (A.9) and carrying out straightforward operator algebra, we 
obtain step by step 

[n,V] = jdki j dk2 [F^R'^'F, FW^^F] a^{k^)a^{k2) 

= jdk, j dk2 F\R^\V^^]Fa,{ki)a,{k2) , (A.12) 
where in accordance with Eq. (A.8') 

F^[R^\V^']F = j dp' j dp Y.F\p',r')[R^\V^'%p'r'-pr) F{p,r) , 

r'r 

and it is implied that 
[Rki^V''^]{p'r';pr) 

^ fdq J2[R^'{p'r'; qs)V'''{qs; pr) - V^^{p'r'- qs)R^'{qs-pr) ] . (A.13) 

■' s 
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In the above calculations it has been convenient to employ the identity 

[AB, CD] = A{B, C}D - {A, C}BD - C{D, A}B + CA{D, B} 

for four operators A, B, C and D. 

Further, after applying another useful relation 

[AB, CD] = A[B, C]D + [A, C]DB + AC[B, D] + C[A, D]B , 

a similar derivation for [7^, V^] results in 

[7^, V^] ^ Jdki J dk2 {Ft[i?'=i, T/-'=2]F4(fc2)ae(fci) + 5{ki - k2)F^R^'F ■ F^V-^'F} . 

(A.14) 

Now, we see that the g'^-order commutator [R, V] brings in the total Hamiltonian K, 
see Eq. (2.23), the interaction terms which describe the following real physical processes: 

a) the TTTT — > // creation and the // tttt annihilation from [7^, V] (see Eq. (A. 12)) and 
its H.c; 

b) the tt/ nf, nf ^ tt/, // ff, ff ff and // ff scatterings from [7^, V^] 
and its H.c. (see (A.14)). 

In addition to these contributions, [R, V] contains interactions, which have nonvanish- 
ing matrix elements between the vacuum Q and two-particle states (e.g., those responsible 
for the virtual process Q — > tttt), and between Q and four-particle states (e.g., for transi- 
tions Q — > tttt//). There are also interactions responsible for the transitions one particle 

three particles and one particle one particle. Except the latter, all these terms are 
"bad", i.e., they hinder Q and one-particle states to be K eigenstates. 

For example, let us consider the term dcdlacttc which enters in [7^, V]. It has nonzero 
matrix element {fl \ dcd^a^ac \ tttt) that becomes evident after normal ordering, 

dc{-p',r')dl{-p,r)ac{ki)ac{k2) = -dl{-p,r)dc{-p',r')ac{ki)ac{k2)+Sr'rS{p'-p) 00(^1)00(^2) • 

This illustration shows that normal ordering is a constructive tool in the framework of 
our clothing procedure. As in many applications of the method of second quantization 
(e.g., in field theories of the evolution operator or the S'-matrix), this operation enables 
us to classify the separate contributions to the original Hamiltonian at every stage of the 
clothing procedure. 

A. 4. Now, we shall discuss in detail the origin of two-operator meson terms which 
stem from the commutator [R, V] . They are essential elements in our treatment of 
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the particle mass renormalization (see Sect. 2.5). These terms appear after reshuf- 
fling the operators of the expressions F^[R''\V''^]Fa,{ki)ac{k2) (see Eq. (A.12)) and 
F''[R^\V^^^]Fal{k2)ac{ki) (see Eq. (A.14)) into normal order. 

In the first case, it touches upon the terms of the dd^-kind and leads to 

F^R^\V^']F =: F^[R^\V^']F : +Tt[R^\ v''']22 , (A.15) 

where 

Tr[R''\V%2 = [dp Y:[R''\V%2{pr;pr) , (A.16) 

[R^\V^^]22{p'r']pr) 

= [dq Y:[R2l{p'r';qs)Vp2'{qs;pr)-V2f{p'r';qs)R%^{qs;pr) ] . (A.17) 

The trace Tr is evaluated using the properties of the 7-matrices and Dirac spinors. We 
find with the solutions (A.l) 



^fc.]^^ = 2-^5(fei + k2) , (A.18) 



_ f. | Sp[A-(-p)A_(p-fc)] Sp[A-(-p)A+(-p + fc)] 

^ Sp[A_(-p)A_(p-fc)] ^ Sp[A_(-p)A+(~p + fc)] | 

with the notations A+(A_) for the projection operators onto the nucleon positive 
(negative)-energy states: 

= ^ ■ 

While deriving Eq. (A. 19), we have taken into account that 

Y,v{ps)Ov{ps) ^ -Sp{A4p)0} , 

s 

where O is a combination of 7-matrices. 

Of course, one can collect the similar terms with the same numerators in the r.h.s. of 
Eq. (A. 19) ( e.g., the first term with the third one inside the curly brackets). However, 
we prefer other continuation that enables us to get immediately an explicitly covariant 
form of tj^. 

First of all, we find 

EpEp_k + p- {p - k)+m'^ 



Sp[A_(-p)A_(p-fc)] 
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Substituting these expressions into (A. 19) and uniting therein the first term with the 
second one and the third term with the fourth one it can be shown that 

^2 



9 



dp 



pk 



2(2-TTf J Ep [i^'^ + 2p_k ii'^-2pk 



or 



9' 



dp 



1 + 



2(27r)3y Ep 

with p_ = {Ep, -p) , p = {Ep,p) and k = {ujk,k). 
Since this integral is a Lorentz-scalar, one can write 

^k — ^fc=o 



(A.20a) 



(A.20b) 



9^ 



2(27r) 



E: 



P 



or 



tk 



47^2 



\/ — /x^ 
/2 arctan 



4£;2 - J 



where the integrals 



and 



X 



^4m'^ — 



(A.21) 



(A.22) 



dx 



dx 



are, respectively, quadratically and logarithmically divergent. 

Now, the resulting contribution to [TZ, V] which is of the OcOc-kind can be written as 

2 / dfci / dfca S{ki + fca)— ac(fei)ac(fc2) = 2 [ dk — ac(fc)ac(-fe) . (A.23) 
J J a;^j J cuh 



^k 



The commutator [R, V] includes also the Hermitian conjugate of [TZ, V] . Therefore, we 
obtain from [JZ, V] + H.c. the following "bad" two-operator meson contribution: 



2 [dk ^{a^{k)a^{-k) + al{k)al{-k)} 



(A.24) 



In the case of F^[R'^\V-'^']Fal{k2)a^{ki) (see the beginning of this subsection) after 
normal ordering of the dd^-kind terms one has to deal with Tr[R^'' ,V~^^]22 that differs 
from Eq. (A. 16) only by the replacement ^2 — — ^2- Therefore, we obtain from [7?., V^] 
the following expression bilinear in the meson operators: 



2 [dk ^al{k)a,{k) . 



(A.25) 
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The same expression stems from the H.c. of [7^, V^], 

Finally, uniting all these results one can write the entire contribution from | 
K{a), which is bilinear in the meson operators, 



[R, V] to 



dfe -^{2al{k)a,{k) + a,{k)a^{-k) + al{k)al{-k)} 



(A.26) 



Appendix B 



Mathematical Aspects of the Clothing UT. 



By way of a simple example we shall show that W used in the framework of "clothing" 
approach (see Sect. 2) may happen to be not an unitary operator in the usual sense: 
unitary operator transforms vectors of a Hilbert space into vectors of the same space, the 
scalar products being conserved. We argue that W need not be such an operator, viz., 
the clothing program can be described using an algebraic language as if W is an element 
of some algebra, being unitary in an algebraic sense. Our example shows that such W can 
have a representation by an operator that transforms vectors of a Hilbert space Tio into 
vectors of another Hilbert space H which is orthogonal to Hq. In general, the operator 
representation of W turns out to be unnecessary because W is not used in calculations 
of probability amplitudes, expectation values and other quantities which have physical 
interpretation. We note that "clothing" allows us to choose a proper Hilbert space for 
field model with the total Hamiltonian H = Hq + Hj. This space usually is different from 
the space spanned on Hq eigenvectors. 

B.l. The clothing program can explicitly be carried out in the following model. The 
scalar (meson) field (j){x) interacts with a fixed (external) source, the Hamiltonian being 

= H^^^ + 9 1 4'{x)f{x)dx, where H^'^^ is the free meson Hamiltonian (see Eq.(2.1)). 
The Hamiltonian can be represented as 




(B.l) 



with 





9 



where o;^ = y fc^ + /Xq is the bare meson mass. The Fourier transform / (k) of the source 
function f{x) is a constant if the source is pointlike: f{x) ~ 6{x). 
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In accordance with the recipe (A. 4) we find the corresponding clothing transformation 
W = exp R = W{ac) = W{a) with its generator 

J v{k)[a^{k) - a{k)]dk (B.2) 

When deriving this expression, it is convenient to exploit the relation 

For brevity, we suppose that the form factor v{k) is real {v*{k) = v{k)). 
Further, one can verify that 

a(k) = Wac(k)W'' = a^k) - v(k) (B.3) 

and 

H = K{ac) = y u;laK{k)ac{k) - L , (B.4) 

where L = J v{kYdk is the c-number that shifts the H spectrum. 

Thus, this transformation W not only removes from H the "bad" terms (hnear in a 
and a)) but it reduces the primary eigenvalue problem to a very simple one. In our model 

K{ac) does not contain interaction terms: clothed particles (mesons) turn out to be free. 
Now, we are interested in the clothed no-particle state, i.e. the vector fl such that 

ac(k)n = 0, Vfe (B.5) 

with 

(fi I fi) = 1. 

One can show that the vector W'^{a)ilo obeys the condition (B.5). In fact, 

ac(k)W\a)no = W^a(k)no = (B.6) 
So, one can put fl — W^{a)flo, i.e., 

n = exp{- J v{k) [a^{k) - a{k)]dk}no = 

exp{-L/2)cxp{- J v{k)a\k)dk}no (B.7) 
Here, we have used the Hausdorff-Weyl formula, 

for the two operators A and B such that [^4, B] commutes with A and B. 
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Eq. (B.7) represents the "clothed" vacuum Q as the superposition of "bare" states 

Qo, at(fci)fio, at(fci)at(fc2)l^o,... (B.9) 

Expansions similar to (B.7) can be written for "clothed" one-particle state a^c{k)fl as 
well as for all vectors of the kind 

W^V(fc,) . . . at(fc„)fio = at,(fci) . . . aKikn)W^no , (B.IO) 

each of them being the H eigenvector. 

The states (B.9) are the Hq eigenvectors and they form the basis of the Fock (Hilbert) 
space Hq. The set 

a^c{ki)n, a^c{ki)a^c{k2)^,... (B.ll) 

of the H eigenvectors is the basis of the Hilbert space which we shall call Ti. We see that 
in our model with the cutoff function v{k) which decreases rapidly enough to yield a finite 
normalization factor exp(— L/2), the space Ti. can be spanned onto the vectors (B.9). 

It is not the case in the model with a "soft" form factor v{k). Indeed, if L = cxd 
(this occurs, e.g., for the pointlike source) then we obtain the zero values for all VT^fio 
projections on the vectors (B.9). Moreover, all vectors (B.IO) are orthogonal to Tio if 
L = oo. We may conclude that all vectors (B.ll) are zero if Tio is assumed to be the 
complete space of states. The operator W"^ transforms all T^o vectors into the zero one 
and therefore can not be called the unitary operator because the latter, by definition, 
must transform vectors from TYq into vectors from TYq and conserve the scalar product in 
the space. 

Note that in the case with L —>■ oo there are different ways for calculating the scalar 
product (fi I ri) = {W'^Qq I W^'Qo). On the one hand, putting L = oo in Eq.(B.7), we 
obtain ly^fig = once Tig is complete. Therefore, {W'^Qq \ W'^Qq) = 0. On the other 
hand, calculating at first {Q \ Q) at finite L, we obtain unity for any L, i.e., such limit of 
{W^no I W^no) as L ^ oo is unity0. 

B.2. So, in the case L = oo the standard approach that is relied upon the initial 

Hilbert space Tig does not allow to find H eigenstates (B.ll). The situation was clarified 

by Van Hove He considered H as an operator in the space constructed as the 

infinite product of the Hilbert spaces for the oscillators ookd^lk) a{k) (Van Hove assumes 

that k runs discrete values). This extended space is not a Hilbertian one, but it can be 

^°In the context, we find in ||l^(Ch. XII) that (Q \ Q) ~ oo at L — oo. It follows from the supposition 
made therein that the projection ^'^^^ — {flo \ ^) is not zero at L — oo. However, this supposition is 
wrong ( see Eq. (B.7) and Van Hove paper p9[) 
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decomposed into a direct sum of mutually orthogonal Hilbert spaces, Tio being one of 



them (e.g., see |]50|-||52 



By using this space Van Hove has proved that the H eigenvector which belongs to 
the least H eigenvalue (it coincides with the no-particle state Q because of Eq.(B.4)) has 
unit norm and is orthogonal to Tio- Moreover, all H eigenvectors (B.ll) are orthogonal 
to Tio- The vectors (B.ll) form the basis of another Hilbert space which is orthogonal to 
TCq. We have called it H above. 

Of course, the model under discussion is too simplified and the resulting theory is 
equivalent to the free one, viz., H contains no interaction terms after "clothing" ( see 
Eq. (B.4)). But it enables us to suspect that in nontrivial cases the total Hamiltonian 
eigenvectors may happen to be orthogonal to the initial Hilbert space Hq, i.e., "bare" 
states space or space of the eigenvectors of the free part Hq of the total Hamiltonian. 
On this ground one may cast doubt about validity of the usual quantum postulate that 
H as well as other observables can be defined initially as operators in Tio- One can also 
anticipate that the "clothing" transformation may not be unitary operator in Tig- 

B.3. In order to avoid the troubles described above we suggest the following algebraic 
approach to the quantum field theory. We consider all the operators occurred in quantum 
field theory as elements of an algebra, devoid of operator representation. This means 
that ab initio we do not introduce the notion of vectors, describing states of the physical 
system. In this algebra besides the addition and multiplication of elements an operation 
of involution f is defined, which corresponds to the Hermitian conjugation in the operator 
language (see, e.g., [Q (§1-5) and (Ch. HI)). The algebra contains A'^ along with the 
element A. All elements of the algebra (in particular, the Hamiltonian) can be expressed 
in terms of some basic algebraic elements. In the case of Yukawa model the latter are 

a{k),a\k),b{p,r),b\p,r),d{p,r),d\p,r) yk,p,r (B.12) 

The multiplication operation is noncommutative: AB may not be equal to BA. The 
commutators [A,B] for basic elements are assumed to be known (see, e.g., Eq.(2.3), then 
the commutators of any two elements can be calculated. 

The "clothing" program can be formulated and realized using this algebraic language. 
For example, the requirements i) and ii) from Sect. 2 can be replaced by equivalent ones: 
the Hamiltonian H when expressed in terms of new basic elements ap (instead of the 
starting elements denoted as in Sect. 2) must not contain "bad" terms (see Subsect. 
2.4). The elements ap and ap are connected by the isomorphic transformation ap = 
W^ttpW, W being a unitary element of the algebra, i.e., such that W^W = W = 1. Let 
us stress that ap W is calculated using purely algebraic means, namely commutation 
relations and Eq.(2.15). 
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Of course, our theory must provide numbers, which can be compared with exper- 
imental data (cross sections, expectation values, etc.). In quantum mechanics this is 
accomplished by means of the realization of the above algebraic elements as operators in 
a vector space. We define the space as follows. 

After "clothing" we introduce the additional notion of a state Q (cyclic state) such 
that 

ac{k)^l = bc{p, r)^l = dc{p, r)n = Vfc, p, r 

This state coincides with the H eigenstate corresponding to the least H eigenvalue. We 
assume that the observable no-particle state is described by Q. By assumption observable 
one-particle states are described by the states 

a^c{k)Vl, b^c{p,r)n, d^c{p,r)n 

The vectors 

a'c{ki) . . .a^^cikn)^, 6'''c(pi, n) . . . 6"^c(p„, r„)fi, . . . n>2 

can be chosen as the remaining basic vectors of our Hilbert space Ti. Using the vectors 
one can calculate quantities of physical interest in the usual manner. 

One may say that one of the goals of the "clothing" is to select a Hilbert space 
which would be suitable for the given field Hamiltonian H. The first-quantized quantum 
mechanics uses only one universal Hilbert space for different Hamiltonians, the second- 
quantized theory needs distinct spaces for different interactions (e.g., for different values 
of the coupling constant g). 

The approach described above needs not initial space TCq. We need not to consider 
W as an operator in a space. Algebraically W is defined as an unitary element of the 
described algebra. Besides W the algebra has other elements which need not operator 
realization, the free part Ho of the total Hamiltonian being the example. 

Appendix C 

The UT Method in Scalar Field Model 

In order to illustrate the key points of clothing and Okubo's approaches let us consider 
another exactly solvable model 0, in which a neutral scalar field is coupled with 
spinless fermions whose energy is independent of momentum. The model Hamiltonian is 
H = Ho + V, 

H = moB{0)+ [ uJi,a^{k)a{k)dk, ojj. = ^ k^ + fi'^ , (C.l) 
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V^g [ujk [B{k) a{k) + H.c] h{k^)dk, /i(fe^) = , -^^^ ^ (C.2) 
J V2(27r^ 



with 



B{k) = jh\p+k) b{p)dp = B^{-k) , 

where a{k) and b{p) are the destruction operators for bosons and fermions, respectively, 
which meet the usual commutation rules (cf. Eqs. (2.5)): 

[a{k),a^{k')]^S{k-k') , (C.3a) 

{b{p),b\p')} = 6{p-p') . (C.3b) 

The translational invariance of the Hamiltonian provides the momentum conservation. 
The cut-off factor f{x) is assumed to fall off rapidly enough for large x to make finite all 
the integrals that occur in the theory. 

C.l Again the clothing transformation of this Hamiltonian can be found in a closed 
form (cf. App. B). Indeed, calculating the respective integral (A. 4) and noticing that 

[B{k),B{k')] =0 V fc , fe' , (C.4) 

we get 

R = -g j [B{k)a{k) - H.c] /^(fc^) dk = g{X^ - X) (C.5) 

with 

X^ J h{k'^)B{k)a{k)dk . 
Further, it is useful to keep in mind the following formulae 

[X, Xt] = J h\k^)B{k)B'^{k)dk (C.6) 

and 

[X,[X,Xt]] = 0, (C.7) 
whence due to Eq. (B.8) the transformation of interest can be written as 

W = exp[^(Xt - X)] = exp(^Xt) exp(-^X) exp(-^[X, X^]) . (C.8) 
Now, by using the relation (cf. Eq. (B.3)) 

a{k) = Wa^{k)W'^ = a^{k) - gh{k^)B\{k) , (C.9) 
where the boson-type operator 

B,{k) = j bK{p + k) b,{p)dp (C.IO) 
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commutes with W = W{a, h) = W{ac, be) (see Eq. (C.4)), one can show that 

H = H(a, h) = X(ae, 6c) = WH{a,, K)^^ 

= mo5c(0) + j uJkaK{k)ae{k)dk - j u;kh\k'^)B^e{k)Be{k)dk . (C.ll) 
Reshuffling the fermion operators in normal order in the r.h.s. of Eq. (C.ll) we obtain 

J ujkh?{k'^)B\{k)Be{k)dk = Se(0) j u}kh\k^)dk 

- jujkh\e)dk JdpJ dp' bKip + k) bK{p') b,{p) b,{p' + k) . 

The first term in the r.h.s. of this equation has the same structure as moSc(O) in Eq. 
(C.ll) and gives the radiative correction (renormalization) to the bare fermion mass mo, 



m 



^rriQ-g^ j ujkh\k^)dk . (C.12) 



Meanwhile, the H eigenvector belongs to the H eigenvalue m. 

So, 

H = K{a,,b,) = KF + Ki, (C.13) 

Kp = mSc(O) + j u;ka^c{k)ac{k)dk = Kf^rm + -f^ftoson , (C.14) 

Ki = Jdx jdx'^p\{x) ij\{x') Vffd a; - a;' I) ij,{x) ij,{x') , (C.15) 

Vff{\ X I) = Ukh\k'')e^^^dk , (C.16) 

where in agreement with the secondary quantization prescriptions we have introduced the 
ipc - field for clothed fermions in the Schrodinger picture assuming 

^^;e{x) = -jL= [ 6e(p)e^P^dp 

{Mx),tljK{x')}^S{x-x'). 

Therefore, Vff{\ x |) can be considered as a two-fermion interaction potential. 

One should point out that the new interaction Hamiltonian Kj expressed through 
clothed operators no longer contains any self-interaction and leads merely to an interaction 
between pairs of fermions. 

The Kp eigenvectors 

fi, ate(fei)fi, b^ciPi)^: at,(fei)ate(fe2)fi, b^ c{Pi)b^ aK{ki)bK{p^)^l, . . . (C.17) 
with the running element 

Oe(fei . . . fe,; . . . pj = at,(fci) . . . aK{kr)bK{p^) . . . bK{p,)Q (r, s = 1, 2, . . .) (C.18) 
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form a basis (see App. B). 

Another basis can be composed of bare vacuum f^o and vectors 

Qo(fei ...K;p,...p,) = at(fci) . . . a^iKpip,) . . . b^ip,)^^ (r, s = 1, 2, . . .) , (C.19) 
i.e., of the Ho eigenstates. 

Note also that in the given model, due to the absence of // -processes, the vacuum Q 
for the coupled fields coincides (to an accuracy of a phase factor) with the vacuum fig for 
the free fields. In the context, with the help of Eq. (C.9) we find, 

a^c{k)^ = a^c{k)^o = a^(fe)l^o , (C.20) 

i.e., the bare and clothed one-boson states are the same. 

As to the clothed one-fermion states b'^dp)^, they are complex superpositions of 
vectors (C.19) with one fermion and arbitrary boson configurations (a boson cloud sur- 
rounding the fermion). This follows directly from the explicit expression for 

fet^(p) ^ W^b^{p)W = J^iP- P')bKp'W , (C.21) 

J^{q) = .^-L_ J e-'<l^exp{-g J exp[-tkx] (at(fe) - at(-A;)) h{k^)dk}dx . 

The factor J^(g) characterizes a boson distribution in the cloud. In the free case with ( 
^ = 0) J'{q) = 6{q). 

C.2 It follows from Eqs. (C.13)-(C.15) that clothed mesons do not interact with 
nucleons, they are free. In particular, no-meson states cannot pass in states with one, 
two, etc. mesons. This means that K{ac, be) possesses the property of the operator 
H' — WHU obtained from H via Okubo's transformation considered in Subsect. 6.5., 
viz., — H21 — 0. In other words, H' has vanishing matrix elements between the no- 
meson states and one-, two-, . . . meson states. This property of K{ac, be) is specific for the 
given simple model. In general, the clothing UT permits K to contain such nondiagonal 
terms that are responsible for the process NN — > NNn. 

Since the model clothing UT W fulfills Okubo's requirement it is interesting to compare 
the block structure of IF^ with that described in Subsects. 6.2 and 6.3 for Okubo's UT. 
To do it in a compact form one may calculate separate blocks of W"^ decomposition of 
the type determined by Eq. (6.14 ). In the context, recall that Okubo's UT can be 
represented by a function U (a, b) of bare creation (destruction) operators. W"^ can also 
be given by a function W'^{a,b) of the bare operators since W{ac,bc) = W{a,b) (cf. Eq. 
(2.13)). Therefore, for Okubo-like decomposition of one may use the projector r]i 
onto bare no-meson states. 
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In order to carry out the comparison with a solution of the corresponding decouphng 
equation (6.22) (see C.3 below ) we shall find the operator defined by 



(C.22) 



It determines the basic element A"^^ of Okubo's UT (see, for instance, Eqs. (6.8) and 
(6.20) ). 
We have 



W^^^ — rjiW^Tji = Tjiexp (^—gX^^ exp I — — 



exp (gX) rji 



exp -- 



X, X^ 



(C.23) 



lyt^i = 7]2W^7]i = 7]2exp (-gX^) exp [x, X^]j rji = T^sexp (-gX^) 7]iW^^^ 



While deriving these formulae we have used the property 

ac{k)7]i = Vi(iK{k) = Vfe 

and its consequence 

ex.p{gX)r]i ^ rji ^ r]iex.p{-g X^) 
We also have taken into account the relation 



X, xt' 



Vi = Vi 



X, x^' 



(C.24) 
(C.25) 
(C.26) 

(C.27) 



It follows from (C.23) and (C.24) that 

A'' = mew {-gX^) m = (exp {-gX^) - l) r/i (C.28) 

Note that we have managed to find the operator A^^ in an explicit form without solving 
Eq. (C.22). 

Unitarity relations for the Okubo-type blocks = {W'^y^{i,j — 1,2) with C/^^ = 

C/llt(l + ^21t^21)^ll ^ ^llt ^ ^21 (^-^ ^ J21^ ^11 ^ ^ (<^_29) 

^ _^21t[;22 ^ (C.30) 

and 

U^^\1 + A'W^)U'^ = V2 (C.31) 
can directly be verified. For instance, we get step by step, 

_^2it^22 = (1 _ exp {-gX)) exp \g{X - X^)] tj^ = 
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r]iexp [g{X - X^)] 7]2 - r/iexp (-gX) exp [g{X - X^)] r/2 = W^^^ = U^^ , (C.32) 
since in accordance with Eq. (C.27) 

77iexp {-gX) exp g{X - X^) t]2 = r/iexp (^-gX^^ exp — X, 



V2 



T^iexp 



V2=0 



and 



r?2exp 



^22t(i + ^21^21t)f;22 _ f;22tf;22 ^ ^12t^ 



12 



r72exp 



g{X-X^)] 7]2 + V2ew[g{X^-X) 
ri2WW^ri2 = r]2 



Tji exp 



(7(X-Xt) 



^?2 = 

(C.33) 



Q.E.D. 

C. 3 Let us try to solve Okubo equation (6.22) for this model Hamiltonian. It is of 
great interest to attain some experience in handling with similar nonlinear equations. We 
prefer to start with the equation equivalent to Eq. (6.22) : 



for the operator Ji = (1 + ^^i)^^ Eq. (29)' in [0]). 
One should point out the properties 



(C.34) 



Ji(l-r/i) = 



and the condition 



(C.35a) 
(C.35b) 

(C.36) 



Further, introducing the interaction picture operators O, 

0{t) = exp{tHot)Oexp{-%Hot) , 

and noticing that rii(t) = rji one can get the Riccati-type differential equation for Ji(t) 
(cf. Eq. (31) in [0 ), 



tV2-^Mt)% = V2{Vit)J,{t) - Ji{t)Vit)J,{t)}r]^ 
at 



(C.37) 



In the scalar model we have 



V{t) =g uJk \B\k) a^{k, t) + H.c.l h{k'^)dk 



(C.38) 
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with 



so that 



i{k,t) = a{k)exp{~tu!f^t) 



[a{k,t),a\k',t')] = 5{k - k')e^^[-iujk{t - t')] 



(C.39) 



Trying to solve Eq. (C.37), one should note the relation 

V{t) = ~^gJ^ [X\t)-X{t) 



(C.40) 



where 

X{t) = J h{k'^)B{k)a{k,t)dk . 

Now, applying the Lagrange method well known in the theory of ordinary differential 
equations (see, e.g., |Q, §1.6) let us search a solution of Eq. (C.37) in the form 



Ji(t) =exp [~gX^{t)\ G{t) 
where in accordance with Eqs.(C.35)-(C.36) we have 

y]iG{t) = r]i , 
G(t)(l-r/i)=0 

Differentiating (C.41) and taking into account that 

|xt(t),Xt(t) =0, 

we find 



(C.41) 

(C.42a) 
(C.42b) 



f^Mt) = -g (l^n^)) exp [-gX^it)] G(t) +exp [-gX^{t)] jG{t) . 

Substitution of this expression into Eq. (C.37) enables us to remove in the r.h.s. of this 
equation not the linear term V{t)Ji(t) as a whole but its part —ig-^X\t)Ji{t) (see Eq. 
(C.40)). So, we get the following equation for the operator function G{t) , 



?72exp 



-gr]2exp 



-gX^{t)] jG{t)v, = gV2exp [-gX^{t)] {jX{t)+g 



-gX\t)\G{t)r^,{jX{t)+g 



(C.43) 



At the point we have used the properties r^iexp —gX'^(t) = rji and r]i-^X'^(t) = and 
the relation 



exp 



gX^it)] ^X(t)exp[-^Xt(t) 



d 
dt 



Xit)+g 



XKth^Xit) 



61 



that follows from Eqs. (2.15) and (C.4). 
Further, we have 



z J (jjkh'^{k'^)B^{k)B{k)dk = C 



In other words, this time independent commutator acts merely on the fermionic degrees 
of freedom. Note that C commutes with the projector rji. 
Eq. (C.43) can be satisfied if we put 

-gG(t){j^X(t) + aC}G(t) . (C.44) 
It may be shown that one of its possible solutions is 

Gl{t)=r|^ (C.45) 

It is evident that Gi{t) meets necessary requirements (C.42) 
The corresponding operator Ji(t) — exp —gX^(tj^ rji yields 

A''{t) = {exp [-gX^it)] - l}rii , (C.46) 

that is equivalent at t = to the result (C.28). 
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Figure 1: Schematic representation of separate contributions to the effective operator 
Ki. a) hc^hc^hchc, b) djdjdcdc, c) bjajbcttc, d) djajdcttc, e) bjdjacac, f) ajajdcbc, g) 
bjdjdcbc 



Pi ^ ^ P2 Pi ^ ^ P2 

a h 



Figure 2: The g^^-order Feynman diagrams for ttN scattering: a - the s-pole graph; b 
the u-po\e graph 
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Figure 3: The one-pion-exchange Feynman diagrams for NN scattering 
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Figure Captions 

Fig. 1. Schematic representation of separate contributions to the effective operator 

a) bjbjbcbc, b) djdjdcdc, c) bjajbcac, d) djajdcac, e) bjdjacttc, f) ajajdcbc, 
g) bjdjdcbc 

Fig. 2. The g'^-oidei Feynman diagrams for ttN scattering: 
a - the s-pole graph; b - the w-pole graph 

Fig. 3. The one-pion-exchange Feynman diagrams for NN scattering 
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